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'All linear systems are equally happy, whereas each nonlinear sys- 
tem is unhappy in its own way. ' 

Yu. Bolotin, A. Tur, and V. Yanovskii 
"Constructive Chaos" 



Preface 



Nonlinear dynamics is one of the most important and prospective trends of the 
development of economic science. Powerful modern techniques of qualitative 
theory of differential equations and related subjects of mathematical topology 
provide broad possibilities of obtaining substantial results of qualitative charac- 
ter, in the first place, in solving the problems of economic forecasting. Mathe- 
matical economics is characterized by two principally different ways of modeling, 
i.e., static and dynamic ones. Following N. D. Kondratev [25], we shall dwell 
on a more detailed and concrete characterization of these two purely theoretical 
approaches to the study of economic reality. 

The static theory considers economic processes in terms of their instant man- 
ifestation, without any regard to inertial changes in time. The static approach 
to the modeling of economic reality is based on the concept of the equilibrium of 
interrelated elements of an economic system. The concept of equilibrium itself 
had been sufficiently well familiar to scientists concerned with mechanics before 
the appearance in 1776 of Adam Smith's prominent work "An Inquiry into the 
Nature and Causes of the Wealth of Nations" wherein, as it seems, the author 
managed to find an analogy between the economic balance and the resultant 
force in mechanics. A. Smith put forward the most substantial postulate of the 
general theory of equilibrium: namely, an ability of the system of competition 
to achieve such a distribution of resources that, in a sense, proves to be efficient. 

On positions analogous to those of A. Smith stand also the constructions of 
D. Ricardo who put forward the thesis of the freedom of competition and the 
freedom of movement of labor and capital from one sphere of economic relations 
to another. D. Ricardo was aware of the fact that the actual price, the actual 
level of wages and incomes are variables with respect to time, and, at that, he 
suggested that there exist a tendency towards attaining a certain natural level 
of equilibrium for the above mentioned characteristics. 

An exhaustive formulation of the general concept of equilibrium rightfully 
belongs to L. Walras, a representative of the school of marginalism (or marginal 
utility). Works by L. Walras, S. Jevons, and V. Pareto unified the theory of 
equilibrium with regard to an application to the spheres of exchange, production, 
capital and money. They were subsequently elaborated by J. Hicks and P. 
Samuelson. On the whole, these works are of rather broad, comprehensive 
character: economy is considered as a set of individual consumers and producers, 
and the number of involved variables is absolutely unlimited. The system of 



vii 



viii 



PREFACE 



general equilibrium is closed in the sense that the whole set of variables is 
determined by given conditions. In order to verify the compatibility of the 
system with the state of equilibrium, one only has to compare the number of 
equations with the number of unknown variables. However, there exists the 
problem of the existence of the point of equilibrium and of its uniqueness. 

The dynamic analysis in economics formed in parallel with economic theory 
itself. A confrontation between the dynamic and static approaches can be easily 
traced throughout the whole history of the economic thought. Apparently, the 
reason lies in principal differences between the understanding of the balance 
of the equilibrium of forces and casual dynamics. There is a vast choice of 
literature references concerned with the theory of economic dynamics. Among 
a large number of different problems of evolutionary economics, the problems 
of economic growth and of business cycles are the most important ones. In the 
treatment of N. D. Kondratev, these distinctions between dynamic processes are 
interpreted as evolutionary (nonrecurrent, irreversible) and wave- like (recurrent, 
reversible). By evolutionary (or irreversible) processes one means those changes 
that in the absence of external perturbative interactions flow in a certain, one 
and the same, direction. Examples of such processes are given by tendencies 
of population growth, increase in the total volume of production, etc. N. D. 
Kondratev terms as wave-like (or reversible) those processes of changes that 
at each point in time have their specific direction and, consequently, change it 
permanently. In these processes, the phenomenon, being at a given moment 
in a given state and changing it afterwards, sooner or later may return to the 
initial state. The processes of changes in the prices of different consumer goods, 
interest rates, the level of unemployment, etc. may serve as examples. 

In our view, more attention should be paid to the issue of slow (low- frequency) 
oscillations in economics, i.e., to the so-called long cycles (waves). N. D. Kon- 
dratev himself, while singling out long waves of economic activities, related 
them to the industrial revolution at the end of the XVIIIth - beginning of the 
XlXth centuries, the construction of railway networks, the dissemination of new 
communication means (telephone, telegraph) and of electric power, as well as 
rapid development of automotive industry [2j. Nowadays, the most wide-spread 
concept is that of five long cycles whose length is approximately equal to 50 
years: 

• the end of the XVIIIth century - the first third of the XlXth century; 

• the second third of the XlXth century - the early 90s of the XlXth century; 

• the end of the XlXth century - the 30s of the XXth century; 

• the 40s of the XXth century - the 70s of the XXth century; 

• from the 80s of the XXth century up to the present time. 

Let us try to understand the very nature of the mechanism of the long cycle. 
To initiate the expansion phase of the cycle, it is necessary to accumulate not 
only inventions but also capital as well as a desire of entrepreneurs to increase 
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investments. For the industrialist, the dynamics of profit is a factor of primary 
importance. As a matter of fact, the expansion wave of the long cycle devel- 
ops as a system of mutually related and mutually stimulating phenomena: an 
innovation provides a possibility to improve the production conditions, to re- 
duce production costs and to increase profit, which stimulates entrepreneurs to 
introduce innovations under the condition of availability of necessary resources. 
Innovations give rise to an increase in profit, which generates additional invest- 
ments, an increase in the volume of demand and a general positive movement 
of the growth rate of business factors. 

However, at a certain moment the dynamics of the process exhibits a return 
point. A technological basis for this is provided by substantial weakening of the 
factors that initiated the expansion phase of the cycle. A cessation of their action 
slows down the growth of profitability and then decreases it, which reduces 
the interest of business structures in further innovations and investments. The 
industry slows down its growth, and negative effects of the economic life, typical 
of long-wave decline, appear. In the process of this decline, an increase in 
the number of new inventions takes place, which creates a prerequisite for the 
completion of the decline and the beginning of a new expansion wave. 

N. D. Kondratev's undoubted merit consists in the fact that he based his 
conclusions on an analysis of long temporal series of prices of the commodity 
output, of interest rates, wages, etc. The above-described dynamics is merely 
a simplified scheme, because the actual changes in the economy are much more 
complicated and diverse. 

The modern status of macroeconomics cannot be understood without an 
evaluation of the contribution of J. M. Keynes. As a matter of fact , he created 
macroeconomics as a science in the 1930's in order to explain the causes of the 
Great Depression that proved to be the most large-scale recession of the XXth 
century and, as such, the most important event in the modern history of business 
cycles '18 . The theory, existing at that time, could not explain why the GDP 
of the USA had fallen by one third from 1929 to 1933 and the unemployment 
level had risen to one fourth of the total work force. The classical theories were 
bases on the assumption that the economy was at competitive equilibrium, 
with the market regulating everything. In particular, high unemployment had 
to give rise to the reduction in wage rates down to the level where the employers 
would agree to employ all those who were willing to work for these wages and 
unemployment would disappear by itself. However, in practice, this was not the 
case. Therefore, Keynes put forward radically new ideas whose essence could 
be reduced to two major postulates. 

Firstly, the economy is not at competitive equilibrium at each separate mo- 
ment; that is, the "invisible hand" of the market does not fulfil its duties. The 
basic reason for this equilibrium is the conservation of fixed prices and wage 
rates for a long time and the absence of adaptation to the current market con- 
ditions. Secondly, the level of the development of the economy is determined by 
the aggregate demand, and the latter, in its turn, depends on some unexplained 
factors that were vaguely termed by Keynes as "the brute of investors". On 
the basis of these two assumptions, the edifice of the theory of macroeconomics 
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was erected. According to its postulate, the volume of the GDP of the country 
is influenced by the scale of expenditure of consumers, investors and the gov- 
ernment on commodities and services. Therefore, business cycles are stipulated 
exactly by oscillations of the demand rather than by resources of the country. 
The first and main formal Keynesian model IS-LM was formulated on the basis 
of the Keynesian theory by J. Hicks in 1937. 

In the course of the next several decades, up to the mid-70s, the discus- 
sion went on mostly in the mainstream of the Keynesian theory. The issue was 
whether the government should at all try to revive the economy in the periods 
of decline and, if it were so, by what means. According to the above-mentioned 
position, the government had to react to the decline by increasing government 
spending. From a point of view of other scientists and specialists, stabilization 
should be achieved by means of control over money supply: this point of view 
gave rise to the development of the doctrine of monetarism according to which 
the main objective of the state is to avoid strong oscillations of the money sup- 
ply. Nevertheless, the positions of both the antagonistic scientific trends agreed 
on the point that the basis driving force of the cycle were oscillations of the 
volume of demand, and, therefore, the main differences between "Keynesians" 
and "monetarists" were almost completely obliterated. 

However, in the mid-1970s the word economy faced a new phenomenon, i.e., 
stagflation, that could not be satisfactorily explained within the framework of 
the Keynesian concept. At that time, there appeared critical works of R. Lucas, 
subsequently the 1995 Nobel Prize winner, who criticized not only the economic 
policy of the authorities but also the whole Keynesian theory of business cycles 
for the disregard of optimum behavior of business agents including the forma- 
tion of rational expectations. He suggested that, in contrast to investors and 
consumers in Keynes's models that followed certain formal rules of behavior, 
business agents made, on the average, correct forecasts of the future state of the 
economy and adhered to the strategy of maximizing their own proflts. All this 
created a demand for some alternative theory of business cycles. 

This niche was occupied by the American economists F. Kydland and E. 
Prescott, who won the 2004 Nobel Prize for Economics. In their seminal paper 
pS] . they proposed a new description of a real business cycle based on the fact 
that flrms maximized the profit and made decisions to invest taking into account 
the expectation of future demand for their product and of the development of 
technologies. 

Kydland and Prescott presented a series of dynamic models and showed 
what kind of behavior of basic economic variables (GDP, investment, and sav- 
ings) was to be expected depending on the effect of technological shocks upon 
labour productivity and changes in external market conditions. The authors 
demonstrated that the results of modelling were in satisfactory agreement with 
the observed regularities. Besides, they drew an important conclusion that a 
considerable part of oscillations of GDP in many countries corresponded to the 
predictions of equilibrium models. In other words, there is no need to introduce 
into these models deviations from market equilibrium in the Keynesian spirit 
and to realize governmental stabilization policy. 
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For fairness, it should be noted that the theory of cycles of F. Kydland and 
E. Prescott by no means explains all the phenomena of actual economic reality: 
it is permanently subjected to constructive criticism by " new Keynesianists" . In 
particular, the most striking example of disagreement between representatives of 
these two schools is an attempt to explain the technological boom of the 1990s. 
One is just left with expectations that in not too remote future a consensus will 
be achieved concerning the actual sources of business cycles. 

In our point of view, the achievement of this goal is impossible without ac- 
cepting the fact that economics, in essence, is a developing system and should 
be constructed within the framework of the theory of developing systems whose 
constructiveness is convincingly proved by the example of chemical kinetics, bi- 
ology and other natural sciences. In this theory, it is shown that in the process 
of proceeding to the goal in the presence of substantially nonlinear feedback 
couplings, there emerges a whole hierarchy of instabilities that leads to the ap- 
pearance of limit cycles, homoclinic structures and to spontaneous formation 
of chaos. As a result of such transformations (bifurcations), several different 
states of business equilibrium may appear (the so-called effect of bistability) . 
The methods of nonlinear mechanics allow us to predict the moment of the oc- 
currence of a chaotic regime in the system under investigation, the number of 
possible states of equilibrium and to determine the character of their stability. 
All this, in its turn, generates a principal general problem of the construction 
of alternative scenarios of complex, irreversibly developing system. It would be 
in order to mention here the statement of G. Malinetskii ,28]: "Indeed, social- 
technological objects are complex hierarchy systems, with various processes in 
them developing at different characteristic time scales. The rate of their insta- 
bility, the limits of their predictability are different as well. In the economic 
system, the horizon of the forecast has fallen sharply: whereas just 15 years ago 
5-year directive or indicative planning was a norm in the world, nowadays this 
is out of question. In the world, there is more and more supply of 'quick money' 
and less and less supply of 'slow money'. However, on the other hand, sta- 
ble development of the society requires slowly changing strategic goals, scales 
of social values and norms, culture and ideology. One needs technique, the- 
ory and formalism that would allow one to analyze possible dynamics of such 
'different-time-scale' systems and to direct their development on this basis." 

One can hardly question the fact that exactly the mathematical technique 
of nonlinear dynamics provides the very tools that allow us to approach closely 
the solution of the problem of " designing the future" , of finding stable and safe 
ways of social and economic evolution. The experience of the application of 
methods and models of nonlinear dynamics has shown that many complex de- 
veloping systems can be satisfactorily described with the help of a small number 
of variables, the order parameters. The determination of the order parameters 
is realized by reduction of the multidimensional system to a subspace of a small 
dimension owing to methods of the theory of bifurcations and of the theory of 
central manifolds. However, exactly this fact predetermines the locality of the 
carried out analysis of dynamic behavior of the studied system. Its applicability 
is admissible only in small neighborhood of the bifurcation point, to solutions 
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of small amplitude. In what follows, wc shall present other periodic solutions of 
small amplitude generated as a result of the Andronov-Hopf bifurcation of limit 
cycles and shall determine the character of their stability. 

In this book, the choice of the discussed models is made more or less arbi- 
trarily. The authors consider the models that are rooted in basic principles of 
traditional economics, neoclassical synthesis and Keynesianism. 



Chapter 1 



Instability and cycles in the 
Walras-Marshall model 

Economics operates such notions as the quantity of goods (productive factors) 
and their price. In every market, there exist groups of sellers and buyers. In 
this chapter, a model of the market for one kind of goods will be considered. 
In the model of a single market the variables, i.e., functions of time, are the 
volumes of bought and sold goods as well as their prices. The basic principle of 
modeling of the market interaction is the formation of balance relations between 
the volumes of the demand and the supply of goods and, accordingly, the prices 
of the demand and the supply. 

The problem of joint action of demand and supply as indicators determining 
quantitative relations between the volume of a commodity and its price in a 
given market is very precisely characterized by A. Marshall [29]: "We could ask 
on equal grounds whether the price is regulated by utility or production costs, 
or whether a sheet of paper is cut by the upper or the lower blade of the scissors. 
Indeed, if one blade is kept motionless and cutting is carried out by the motion 
of the other blade, we can, without a good deal of thinking, argue that cutting is 
done by the second blade. However, such an argument is not completely exact, 
and it may be justified only by a pretension to mere popularity rather than to 
an exact scientific description of the realized process." 

For more concrete understanding of the modern phenomenological basis of 
demand and supply, we should present a definition of these notions, using the 
formulations given, e.g., in \W\ . 

By the commodity demand one means the quantity of this commodity that 
an individual, a group of individuals or the population on the whole are ready to 
buy per unit time under certain conditions. A list of these conditions includes 
the tastes and the preferences of the buyers, the price of this commodity, the 
income rate, etc. By the demand price one means the maximum price the 
buyers agree to pay for a fixed quantity of a given commodity. At the same time, 
the dependence of the volume of the demand on its determining factors is called 



1 



2CHAPTER 1. INSTABILITY AND CYCLES IN THE WALRAS-MARSHALL MODEL 



the demand function. 

Analogously, the supply serves as a characteristic of the readiness of the 
seller to sell a certain quantity of the commodity in a fixed period of time. 

By the volume of the supply one means the quantity of a certain com- 
modity that one seller or a group of sellers are willing to sell in the market per 
unit time under certain conditions. 

These conditions, as a nile. include the properties of the applied manufactur- 
ing technology, the price of the given commodity, the price rates of the employed 
resources, tax rates, subventions, etc. The supply price is the minimum price 
at which the seller agrees to sell a certain quantity of a given commodity. The 
dependence of the volume of the supply on the structure of its determining fac- 
tors is called the supply function. Let us point out that the supply function 
as well as the demand function can be represented in three ways: in the form 
of numerical tables, graphically, and analytically. In what follows, we shall use 
only analytical representations for the functions of the demand and the supply. 



1.1 Nonlinearity in the Walras model 

In classical economic theory, one employs two equally admissible but principally 
different versions of the description of the mechanism of an interaction between 
the demand and the supply. The first approach, worked out by L. Walras, 
postulates that the driving force of changes in the price is the volume of excess 
demand under a given instant value of the price. In a dynamic aspect, the 
process of finding the equilibrium in L. Walras's spirit can be represented in the 
form of the differential equation 

^=m(F^(P)-F^(P)), (1.1) 

where P = P {t) is the price of the commodity; 
Y^ = Y^ (P) is the volume of the demand; 
yS _ yS j-p-j jg ^Yie volume of the supply; 
m > is a constant of the time of the limit process; 
t is time. 

The sign of the quantity AY = Y^—Y^, called the volume of excess demand, 
determines the direction of the changes in the price. It is obvious that for 

Ay > the market price rises, whereas for AY < it falls. The condition of 
the existence of an equilibrium price Pe is the existence of the solution to the 
equation 

{Pe) - Y^ (Pe) = 0. (1.2) 

There exists also a different approach to the problem under consideration, 
attributed to A. Marshall. Its essence is that a change of the volume of the 
mass of commodities in a given market is determined by the influence of the 
difference between the demand price and the supply price, to which the sellers 
(or the manufacturers) respond by an increase or a decrease in the volume of the 



1.1. NONLINEARITY IN THE WALRAS MODEL 



3 



supply of the commodity. In a mathematical form, this statement is expressed 
by means of the following differential equation: 

— =m{pD(Y)-P^{Y)), (1.3) 

where Y — Y (t) is the volume of the commodity; 
(Y) is the demand price; 
(Y) is the supply price; 
m > is a time constant. 

In Eq. (|1.3|) , a surplus of the demand price over the supply price stimulates 
an increase in Y; and if the supply price is higher than the demand price, the 
value of Y decreases. An equilibrium value of the volume of the commodity Ye 
is determined from the equation 

(r) - (F) = 0. (1.4) 



The algebraic equations l|1.2p and (|1.3p may have only one or several so- 
lutions. It means that both a unique state of equilibrium as well as a set of 
equilibrium states is possible. It is obvious that the nonuniqueness of equilib- 
rium values of the volume and of the price of the commodity is explained by 
the presence of nonlinear relations in the basic equations. 

An important problem is an analysis of the stability of the available states 
of equilibrium. It is necessary to ascertain the reasons why an equilibrium vol- 
ume of the market remains constant under certain, remaining within certain 
limit values, fluctuations of the price, or on the other hand, why, under a given 
equilibrium price rate, changes in the volume of the commodity also take place. 
In what follows, by the stability of equilibrium we understand an ability of the 
overbalanced market to return again to the initial state owing to the action 
of endogenous factors. Besides, the problem of the stability of market equilib- 
rium is directly related with the problem of the necessity to employ additional 
measures to regulate market relations. 

First of all, let us consider the problem of stability of the economic model 
(|l.ip described L. Walras's theory. Let us set the coefficient m = 1 in Eq. 
(jl.ip . In the neighborhood of the equilibrium point P = Pe determined by 
the solution of Eq. (|1.2p , we can approximately represent the functions of the 
demand Y^ (P) and of the supply Y^ (P) in the form of polynomials obtained 
by the truncation of the corresponding Taylor series 

k k 

Y''{P)^Y.'^{P~Pe)\ Y^P)^Y.^^P~PEr, (1.5) 

dP' dP' ' 

If we introduce a new variable x = P — Pe, which is a deviation of the price 
from its equilibrium value, equation (|l.ip takes the form 

x^Fkix), (1.6) 
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where x = Fk{x) = X;Lo ff^'' = " ^^^^ lEH* follows 

oo = 0, because do — So = Ye , which is an equilibrium volume of the market. 
At the same time, a; = is a stationary point (the state of equilibrium) of Eq. 



The differential equation p.6p is called a dynamic system of the first order. 
The phase space of the considered system is one-dimensional, therefore, the 
studied process of change in the price can be represented by the motion of an 
image point on the phase straight |8j . 

Indeed, in general, the main elements that determine the partition of the 
phase straight into trajectories are the states of equilibrium of the system. The 
values X = xj that make the function Fk (x) vanish are themselves independent 
phase trajectories. The rest of the trajectories consist of line segments between 
the roots of the equation Fk (x) = 0, or of rays forming half- intervals between 
one of the roots and infinity. The direction of the motion of the image point 
along these trajectories is determined by the sign of the function Fk (x): for 
Fk (x) > the image point moves to the right, whereas for Fk (x) < it moves 
to the left. If the form of the curve z — Fk {x) is known, it is not difficult to 
establish concrete partition of the phase straight into trajectories. 

An example of such partition is given in Fig. 1.1, where the arrows show 
the direction of the motion of the image point. From the structure of the par- 
tition of the phase straight into trajectories, it follows directly that the states 
of equilibrium of the system at the points xi, X4 are stable, whereas they are 
unstable at the points X2, X3, x^. It is directly seen in Fig. 1.1 that in the stable 
states of equilibrium the derivative Fj. {x) < 0, whereas in the unstable states 

(x) > 0. The value Fj. (x) = may occur at points of both the stable and 
unstable state of equilibrium. (This situation itself deserves independent con- 
sideration, because it requires some additional conditions for the determination 
of the type of stability of the stationary point.) 

As the character of the change of the variable in the first-order system ()1.6|) 
is completely determined by the explicit form of the function Fk (x), it is of 
interest to consider cases of different values of the order of the polynomial k. 

Let k = 1. Then Fi (x) = aix is a linear function, and there exists the single 
state of equilibrium xe — 0- The stability condition in this case is F{ (0) < 0. 
This inequality reduces to the relation di < Si which is the classical condition 
of stability of L. Walras. 

Let us try to interpret the linear stability of L. Walras using the notion of 
the elasticity of the demand and supply functions to price. 

According to the definition of elasticity, in our notation, we have: 



diPE 



SiPe 



Vd = 



So 



or, taking into account that do — So — Ye, 



diPE 



SiPe 



(1.7) 



Vd = 



Ye 



Ye 
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Figure 1.1: The dependence of the excess-demand function on price deviations. 

where r]jj, rjg are coefficients of the demand and supply elasticities to price. 
They are dimensionless, i.e., relative, quantities. 

Therefore, the inequality ai = di — Si < Q can be easily reduced to the form 

ai - ^ (^D - '/s) < 0. (1.8) 

Given that Yg, Pe are always positive, the condition of L. Walras is for- 
mulated as follows: for the stability of the linear system (|1.6p with A: = 1, it is 
necessary that the elasticity of the volume of the supply to price should exceed 
the corresponding demand elasticity, i.e., rjg > rjjj. In other words, if we intro- 
duce the quantity i] = rj^ ~Vs conditionally termed an excess-demand elasticity 
to price, the stability of (|1.6p is determined by the sign of 77: for 77 < 0, we have 
stability, and, on the contrary, for 77 > 0, we have instability. 

Let us consider the peculiarities of the behavior of the system (|1.6p in the 
case k = 2. Here, F2 (x) — aix + ^x"^ is a quadratic function of the initial 
variable. 

The equation F2 (x) ~ 0, or aix + ^x"^ — 0, has two roots: x\, = Q and 
x\ — — To determine the character of stability of each singular point, it is 
necessary to evaluate F2 {xe)- 

As a result of differentiation, we have: 

F2{xe) ^ ai + a2XE- (1-9) 
The substitution of the values of x\< and x\ in expression (|1.9p yields 
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Figure 1.2: The diagram of the transcritical bifurcation. 



It is thus obvious that the stabihty of both the states of equihbrium is 
completely characterized by the sign of the quantity oi (or 77). 

In this case, if ai > 0, i.e., if the demand is more elastic than the supply, 
the state of equilibrium = is unstable, whereas = "1^7" stable. 

On the contrary, for ai < (the demand is less elastic than the supply), x^ 
is a stable state of equilibrium, and, accordingly, x'^ is an unstable one. 

In a noncoarse situation, when ai is a small quantity changing its sign in the 
neighborhood of zero, the so-called transcritical bifurcation appears illustrating 
a change of stability of the states of equilibrium: see Fig. 1.2, where 02 < 0. 

As a result of this bifurcation, the singular points a;^ and x'^ merge for 
ai = 0, i.e., when the elasticity of the demand is equal to that of the supply, 
to form the single two-fold state of equilibrium xe — 0. At the same time, the 
condition 02 ^ is important. 

Here, we can observe a considerable difference between the behavior of the 
nonlinear system from that of the linear model, which manifests itself in the 
present of two equilibrium values that are transformed into one and the same 
point of equilibrium as a result of a transcritical bifurcation. 

It seems to be reasonable to attribute pithy economic meaning to the co- 
efficient of the quadratic term, 02, in terms of elasticities of the demand and 
supply functions. To this end, it is necessary to find the derivatives of the corre- 
sponding types of elasticity with respect to the price at the point Pe- Skipping 
over intermediate transformations, we present the following expressions for the 
quantities d2 and S2 as functions of 77'^, rj'g, rjjy, rj^: 



Ye 



Vd + Vd) > 




E 




E 



(1.10) 
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Subtracting the second equation of (|1.10p from the first one, we get 
«2 = ^ Hv'd - V's) Pe + {Vd + 7?d - 1) {Vd - Vs)) : 

or 

a2^^{ri'PE + {riD+riD-^)ri)- (I-H) 

^E 

It is worth noting that the dependence of the coefficient 02 on the excess- 
demand elasticity r/ and on its derivative with respect to the price ry' is a linear 
function. 

Let us consider the case when the function of an excess demand in the system 
(|1.6p is cubic. This takes place for k = 3, and, accordingly, 

F3 (x) = a,x + ^x^ + ^x^ (03 ^ 0) . 

The cubic equation (x) = may have, depending on the coefficients, one 
or three real roots, and, accordingly, the system (jl.6p may have one or three 
states of equilibrium. 

By analogy with the previous case, the stability of each state of equilibrium 
is determined by the sign of (a;^). 

Let the system (|1.6p have the representation 

i = aix+^x^ + ^x\ (1.12) 

Assuming the parameters oi, 02 to be small, sign-alternating quantities, 
we consider the deformation of the saddle-node bifurcation with an additional 
degeneracy in the quadratic term [4j. In the saddle- node case, the truncated 
system with ai = 0, 02 = takes the form 

i^^x\ (1.13) 


After some thinking, we can arrive at the conclusion that there exist small 
perturbations of the function ^x^ when the system possesses one or three 
hyperbolic fixed points (the states of equilibrium) in the neighborhood of xe — 
0, as well as certain "unusual" perturbations when the system possesses two 
fixed points, with one of them being nonhyperbolic. All the above-mentioned 
possibilities can be accounted for by adding low-order terms. It is customary to 
represent this deformation in the form of the equation 

u = ei + e2u + ^u\ (1.14) 

6 

where x + ^, 61 = ^ - ai] , 6*2 = ai - t^. 

Here, di, 62 are also small quantities. The dynamics of this vector field is 
formed, with the accuracy of topological equivalence, by its fixed points and the 
types of their stability. Generally speaking, mathematical theory of singularities 
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provides the means for a systematic study of zeros of families of mappings, with 
one of the examples being given by the right-hand side of Eq. (|1.14p . 

Let us find the bifurcation set of parameters in the parameter plane 9i, 02 
by demanding that the right-hand side of ()1.14|) and its derivative with respect 
to the variable u vanish: 

ip{u) = ei + e2u + ^u^ = 0, 

6 

ij' (u) = 02 + Y^"^ = 0- (1.15) 

By eliminating the variable u from both the equations of the system (jl.lSp . 
we arrive at the following bifurcation set: 

Sel + QasOl^O. (1.16) 

The bifurcation diagram of the system p.l4p is presented in Fig. 1.3 for 
the case 03 < 0. The system (|1.14p may possess one or three coarse states of 
equilibrium. These states of equilibrium merge in pairs on bifurcation lines Gi 
and G2 [formula (|1.16p ] that form Neile's semicubical parabola with the origin 
at the point A = (0,0). The point A corresponds to merging of all the three 
states of equilibrium into one state. For the values of the parameters in the 
plane 6*1, 62 that lie inside the parabola (region 2), the system (|1.14p possesses 
three states of equilibrium (two stable states and one unstable state in between), 
whereas for the "outside" values of the parameters (region 1) there exists one 
(stable) state of equilibrium. 

It should be noted that, for bifurcation values of the parameters, a projective 
mapping of the manifold ip {u, 91,62) = onto the parameter space has fold- 
type singularities. The dynamic system in the neighborhood of the bifurcation 
exhibits hysteresis. Let us change the parameters in order to cross the semicu- 
bical parabola, while keeping watch on the stable equilibrium regime (Fig. 1.3). 
In the process of motion from left to right, the "breakdown of equilibrium" of 
equilibrium occurs on the right branch G2, whereas in the case of the reverse 
motion it occurs on the left branch Gi . This phenomenon is called " a hysteresis 
loop". 

From formula (|1.16p . it is not difficult to obtain an explicit form of the 
bifurcation set in terms of the initial parameters oi, 02. Using (|1.14p and (|1.16p . 
we derive the relation 

8aia3 = 3a^. (1.17) 

Furthermore, if we take into account the expressions for ai and a2 in terms 
of the demand and supply elasticities to price, using (jl.Sp and (|l.lip , we obtain 
the following expression for the bifurcation line: 

8P|a377 - 3Ye {v'Pe + {Vd + % " 1) vf ■ (1-18) 

In the equality (jl.lSp . the value of the excess-demand elasticity rj and its 
derivative with respect to the price rj' are small quantities. Therefore, we can 
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Figure 1.3: Bifurcation diagrams of the "triple-equilibrium" bifurcation. 
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argue that the above-mentioned bifurcation can be observed in the initial dy- 
namic system for very close values of the demand and supply elasticities, and 
of their derivatives, i.e., 

Vd ~ v'd ~ V's- (1-19) 

Such a type of behavior cannot be explained by means of the methods of 
comparative statistics traditionally applied in economic analysis. Exactly the 
analysis of the dynamics of the system in the neighborhood of the state of 
equilibrium has shown that the behavior of the system is no longer characterized 
by a unique and smooth reaction to small shifts of the parameters. At the same 
time, there emerge a variety of states of equilibrium, including multiple and 
sudden jumps stipulated by the irreversibility of the flowing processes. 

Analogously, one can carry out an investigation into the dynamics of the 
process of establishing the equilibrium value of production according to the 
concept of A. Marshall, described by the differential equation ()1.3I ) However, 
in contrast to the Walras model, in this case an analysis of the stability of the 
studied dynamic system involves such substantial economic characteristics as 
the values of elasticities of the prices of the demand and of the supply to the 
production volume. 

It should be noted that the above results ask for more profound understand- 
ing of P. Samuelson's principle of correspondence whose validity relies on the 
suggestion of predetermined stability of the economic system, with changeability 
bearing a smooth character. 

1.2 A modified Walras-Marshall model 

Up to now, we have considered the processes of changes of the market price and 
of the volume of commodity production as independent, and the mathematical 
models (|1.8p and (|l.lip have been studied separately. Therefore, in what follows, 
in order to study the dynamics of a model of a certain industrial object, we shall 
make an attempt to unify the equations of Walras and of Marshall into a single 
economic system, where the processes of production and of formation of the 
price are mutually related [13j . 

The issue of price formation in productive economic systems has always 
been and still remains relevant both for theoretical economic analysis and for 
the solution of concrete practical problems of the enterprise as well. In our view, 
it is important to synthesize two major factors of price formation. Thus, on the 
one hand, classical theory of market price formation argues that the market price 
corresponds to the equality of the demand and of the supply in the commodity 
market. On the other hand, according to the theory of the firm, it is known 
that, in the case of production balance, the price of the products manufactured 
by the enterprise corresponds to the marginal production costs. Thus, the first 
approach treats price formation from the point of view of the consumer, whereas 
the second one does it from the point of view of the producer. However, realities 
of the economic life witness that the processes of changes in the price and in 
the volume of production fiow simultaneously and are interrelated. Therefore, 
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it seems to be reasonable to consider the market mechanism of balancing the 
demand with the supply and the production process of accounting for the profit 
and costs simultaneously, within the framework of a unified dynamic system, 
according to the methodology presented in the work [5] . 

As the starting point, we consider a mathematical model describing the 
dynamics of the interaction between the prices and the volume of production 
(manufacturing) : 

aP^D (P) - Y, 

fiY^P^Ps {Y) , (1.20) 

where P is the price of the produced and sold product; 

Y is the volume of the product in natural terms (the market supply of the 
commodity) ; 

D (P) is the market demand for the product in natural terms; 
Ps (Y) is the supply price, equal to the marginal production costs, i.e., 
Ps{Y)^C' (y); 

C (Y) are the production costs (expenses); 

a, P are constant positive parameters describing characteristic times of tran- 
sient processes. 

The first equation of the system of the two ordinary differential equations 
(|1.20[) is the classical model of market price formation in the form of L. Walras 
(or P. Samuelson). They are based on the scheme of price formation searching for 
the balance between the demand and the supply: for D (P) > Y the price rises, 
whereas for the opposite sign of the inequality it falls. The second equation 
of (|1.20p describes the process of establishing the balance between the price 
and the marginal production costs with respect to the production (the value of 
production). Here, it is assumed that the balance is disturbed and it is necessary 
to regulate the volume of production: ii P > Ps (Y), the profit of the producer 
P ■ Y — C (Y) rises with an increasing the production volume, whereas in the 
opposite case one should decrease production activities. 

The model is based on substantial simplifications. Firstly, the production 
is assumed to be single-product. Secondly, a local outlet without competition 
is considered, when the whole supply is formed by one producer. However, in 
spite of the above-mentioned assumptions, the model ()1.20p admits complicated 
types of behavior, and their analysis will be the subject of further consideration. 

A formal analysis of qualitative properties of system (|1.20[) should start with 
the consideration of singular solutions that characterize the states of equilibrium 
of the economic model. 

By making the left-hand sides of (|1.20p vanish, we obtain two constraint 
relations between equilibrium values of the price P* and of the volume Y* : 

D{P*) = Y*, 

P*=Ps{Y*). (1.21) 

Let us assume that the system of algebraic equations (11.201) has, at least, 
one positive solution {P* , Y*). concerning the volume of the demand D (P), we 
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point out that the dependence on the price is a substantially nonlinear function 
and there exists a Taylor expansion up to the third order in the neighborhood 
of the point P* : 

D (P) = do + diiP- P*) + d^^^—^^ + d^^^—^ + O (|P - P*|') , 

where di — '^^J'pi , « = 0, 3. 

The cost (expenses) function is represented by a quadratic function of the 
volume of production: 

C{Y) = Si^ + S^Y + Co. 

where 5i, So, Co are constant parameters. 

Accordingly, the marginal costs (the supply price) are described by the for- 
mula 

C (Y) = Ps (F) = SiY + So. 
Then, the system of equations (|1.21[) can be represented as follows: 

SiD{P*) + So-P* ^0, 

Y* = ^* (1.22) 
Si 

Having preliminarily changed the time scale, it is convenient to study the 
system (jl.20l) in terms of new variables that represent deviations from the equi- 
librium values P = P — P* , Y = Y — Y* . In this case, the system (|1.20l) reduces 
to the form 

p2 p3 
P = diP + d2—+d3—-Y, 

Y^j(^P~SiY), (1.23) 

where 7 = ^• 

As is obvious, equations (I1.23P possess the trivial state of equilibrium _P = 0, 
f = 0. 

To study the stability of the trivial state of equilibrium, we write down an 
explicit expression for the characteristic equation of the linear part of the system 

X' + {jSi-di)X + j{l-diSi)=Q. (1.24) 
The quadratic equation (|1.24[) has negative real parts under the conditions 

jSi < di, 

Sidi < 1. (1.25) 

The inequality (|1.25[) determines restrictions on the parameters of the initial 
system for stability in the linear approximation. 
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Let us consider in greater detail the situation in the vicinity of the boundary 
of the region of stabihty, taking into account the equahty 



where /i is small, sign-alternating quantity. It is obvious that, in this case, the 
divergence of the vector field of the system (|1.23[) is equal to the small parameter 
fi. Therefore, for fj, < the type of the singular point (the state of equilibrium) 
is a stable focus, whereas for /i > it is an unstable focus. In other words, for 
/i = 0, in the neighborhood of the equilibrium state, there occurs the formation 
(annihilation) of a limit cycle as a result of the Hopf bifurcation. 

Let us verify the validity of the conditions of Hopf 's theorem as applied to 
the system (|1.23l) . The eigenvalues are determined (for /i = 0) by the equality 



where = — 1, = j ~ di, i.e., they are purely imaginary. Upon the substi- 
tution of (|1.26p in the quadratic equation (|1.24p and subsequent differentiation 
with respect to the parameter /i, we obtain for /i = 0: 



From (|1.28p . it follows that the real part of the eigenvalue with respect to 
the parameter does not vanish, i.e., the eigenvalues in the complex plane cross 
the imaginary axis with a nonzero velocity. As a result, all the conditions of the 
Hopf bifurcation theorem are fulfilled. 

Let us turn once again to (|1.26p in order to give meaningful interpretation 
of this equality. 

The above-mentioned condition can be fulfilled is the parameters di and 5*1 
have the same sign. As in the following the value Si = C" {Y) will figure as 
the bifurcation parameter, its positivity characterizes concavity of the function 
of expenses C (y), whereas its negativity, accordingly, characterizes convexity. 
From an economic point of view, 5i < determines a positive effect of the 
volume of production (C" {Y) < 0), whereas 5*1 > means that a rise in the 
costs outruns the output of the products (C" (Y) > 0), i.e., the production is 
resource-consuming. 

In order to determine essential parameters of the limit cycle that character- 
ize its stability and the structure of periodic solutions, we reduce the system 
of differential equations (|1.23p to the Poincare normal form by means of a cor- 
responding change of variables P = xi, Y — dixi + ljx2. As a result of the 
reduction, we obtain for /it = 0: 



7C1 = di - ^, 



(1.26) 



(1.27) 



^^A'(0) = l-z^. 
dii ^ ' 2 2uj 



(1.28) 



Xi = —UJX2 + 




dsxl 
6 



did2 x\ didy, x\ 



(1.29) 



X2 = tOXi — 



Lj 2 to 6 
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Using the explicit form of the coefRcients of the nonhnear terms of the system 
(jl.29p . we derive an expression for the first Lyapunov quantity: 



h (0) 



d3 (7 - dl) + 



d\di 



16(7-^2) 



(1.30) 



For li (0) < 0, a stable limit cycle takes place, and a corresponding regime 
of self-oscillations is called "soft". On the contrary, if li (0) > 0, the limit cycle 
is unstable, the self-oscillations break down "rigidly", with a manifestation of 
irreversibility (hysteresis). The case li (0) = is the most complicated one in 
the sense of a variety of phase-plane structures of the system (|1.29|) . because 
there appears a possibility of simultaneous coexistence of two limit cycles (with 
one being stable and the other one unstable) that subsequently merge into a 
single multiple limit cycle. This bifurcation has codimension two and will not 
be studied in detail in this Chapter. 

The periodic solution of small amplitude e (up to a choice of the initial 
phase) is itself written down in the form ^37j 

P {t) = P* +xi (t) , Y{t) = Y*+ dixi (t) + UJX2 (t) , 



XI (t) 



e cos 



27ri 



+ 



e^d2 
12^2 



3di — di cos 



2ijj sin 



Oie' 



X2 [t) = e sm ( — 



+ 



fd2 

12a;2 



idi 



Ant 



2diUJ sin 



. /47rA 

m 



\ T J 



(1.31) 



2 _ 27(Si 



is the amplitude; T (s) = ^ {l + tsE^ + O (e^)), T2 
is the period of the cycle depending, generally speaking, on 



Here, e 

-j^ (8 + 53- 
the amplitude. 

Thus, by the example of the system of two nonlinear differential equations 
(|1.29p . it is easy to establish that, in contrast to a linear system, periodic solu- 
tions are no longer harmonic, and the period and the amplitude of the oscilla- 
tions are interrelated. 

As an illustration of the obtained results, consider examples of economic- 
capacity cycles for different groups of commodities with regard to the depen- 
dence of the demand functions on the income, following the classification of the 
Swedish economist L. Tornquest [361 . 

Example 1. The demand function for essential commodities has a repre- 
sentation 
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which reflects the fact that an increase in demand for these essential commodities 
gradually slows down with an increase in the income and has a has a limit qi > 0. 
The parameter Ci is called the constant of half- saturation of the income. 
Assuming that an equilibrium income is a function of the price, 

E = PD (P) , 

we express the demand in the form D ^ D {P). After corresponding transfor- 
mations, we get 

D(P) = ^1^. (1.32) 

With the help of (|1.22p and p.32p . we write down the equation for an equi- 
librium price: 

{P*f ~ {So - SiCi) P* - Sm - 0. (1.33) 

Differentiating (|1.32p with respect to the price, we obtain the coefficients of 
the demand function: 

^ 51 ^ 6qi , . 

«i = 9, "1 = Q, «i = J- (1-34) 

[P*f {P*f {P*f 

From (|1.34p . it follows that the quantities di, d2 are negative numbers, and, 
therefore, by (|1.30p . the emerging limit cycle is stable. At the same time, the 
reason for the emergence of the cycle is the fact of positive influence of the effect 
of the scale of production of the given group of commodities, i.e., Si < 0. 

Example 2. The demand- for- luxury-goods function is represented in the 
form 

^_ qiD{D-b2) 
D + C2 

By analogy with Example 1, we express the demand as a function of the 
price: 

DiP)^^^^±^, (1.35) 
q2- P 

where 52, ^2, C2 are positive parameters. 

The equation for an equilibrium price has the form 

{P*f - {q2 + So- S1C2) P* - q2 {Sid + So) = 0. (1.36) 

Assuming that (|1.30p has at least one positive root P* , we evaluate the 
coefficients of the powers of P: 

, 92 ib2 + C2) , 2q2 (62 + C2) , 6g2 ib2 + C2) ,^ 

(92 - P*f (<Z2 - P*f (<Z2 - P*)' 

In this case, the coefficients di and ^2 are positive numbers, and the sub- 
stitution of their values in (|1.30p ensures the condition li (0) > 0, which is an 
indication of a catastrophic loss of stability of the limit cycle. As di > 0, the 
emergence of a limit cycle requires the fulfillment of the condition Si > 0, which 
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is possible only for resource-consuming production with an outrunning increase 
in the costs. 

Let us consider one more version of the model (|1.20[) . assuming, as a pre- 
liminary, that the supply price Ps is a nonlinear function of the volume of 
production. For simplicity, we consider the quantities P and Y to be deviations 
from certain equilibrium values P* and Y*. 

Concerning the demand function D (P) and the price, we put forward an 
assumption that they quadratically depend on their arguments, i.e., 

p2 

D {P) = d2— - d,P, 

Ps (Y) ^ - S,Y. 

Then, the system (11.201) can be represented in the form 

P^d2—-diP-Y, 
Y ^ (^P + SiY - S2^y (1.38) 

where 6^ = 7. 

One can argue that a linear analysis of the stability of ()1.38|) completely 
corresponds to the previously obtained results for the system (|1.23p , including a 
verification of the validity of Hopf 's theorem, up to the substitution of for the 
parameter 7. Therefore, assuming the closeness of the bifurcation parameter 
Si to the quantity p-, we transform (|1.38p to the normal form of the given 
bifurcation with the help of the change of variables xi = P, X2 = 
(w^ = b'^ — di). After some transformations and the introduction of a new time 
scale T = Lot, we obtain: 



Xi = —X2 + 



2 



,2 



di {d2~b'^diS2) xj , b^di52 ^2 Qo^ 

X2^xi'i ^ 2 2~^i^2 - J2-;7- (1-39) 

Let us represent the system of two ordinary differential equations (|1.39p 
in the form of a differential equation for the complex variable Z = xi + 1x2, 
Z = xi — ix2'- 

Z^ - Z^ , , 

Z = tZ + g2a-^ + giiZZ + go2-^, (1-40) 

where 511 = (^2 + ^ (^1^2 - 6^ ^2)); 

520 = i^{d2 + 2b^diS2 + i {did2 + b^S2 {b^ ~ 2dl))); 
902 = (^2 - 262^1 52 + ^ (did2 + b^S2 (6^ - 2dl))). 

Now, we possess all the necessary information for the evaluation of the first 
Lyapunov quantity 

'1 = - 2 ■920511- 
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Taking into account the explicit form of the coefficients of Eq. (|1.40p . we 

get: 

- ifc - ■ ^'-^'^ 

From (|1.4ip . it is obvious that for (i| > 6^5*1 the hmit cycle is stable, whereas 
for ^2 > b^S2 the instability of the limit cycle takes place. 

In Fig. 1.4, we present cyclic changes of the price and of the volume of the 
commodity of the considered state of equilibrium for the following values of the 
parameters of the system (|1.38l) : 

b=l; di=0.5; da = 2; S*! = 0.5; S'2 = 1.5. 

The condition = ±6'^S'2 makes the expression for the first Lyapunov quan- 
tity vanish. This may mean that the initial system (|1.38p possesses two limit 
cycles that can merge into one two- fold cycle by means of trajectory compaction. 
This situation is possible, if the second Lyapunov quantity, I2, does not vanish. 
Setting for definiteness 6*2 = as a result of evaluation, we obtain I2 — 0. 
Moreover, ^3, the next (third) Lyapunov quantity also vanishes under the given 
conditions. The fact that all the first three Lyapunov quantities vanish means 
that the state of equilibrium is a center, not a focus. In other words, the initial 
system ()1.38p . for 

S, = ^, S, = ^, (1.42) 

turns into a conservative one, with the conservation of the phase volume. 

Let us write down the explicit form of (|1.38p . eliminating the parameters 
Si,S2 with the help of Eqs. (fO^ : 

P = d2— - diP - Y, 

Y = b^P + diY-^^. (1.43) 

It is important for us to find the first integral of the system (ll.43p . To this 
end, we transform the variables P and Y with the help of the linear substitution 

Y 1 / y\ 1 
n^p-p y..-(p + -)+-, (1,44) 

''l'»e£= = S;,a=-g^.^2 = t'-<if. 

Then, equations (|1.43p can be represented in the form 

n = aY,^ + afYi - {e + 1)^2 + 



F2 = 2aYiY2. 



(1.45) 
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Figure 1.4: The limit cycle in the system (1.38). 
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The system (|1.45p can be reduced to a total differential equation: 

(^aY^ + aS_^Yi - {f + 1)^2 + c?^2 = 2aYiY2dYi. (1.46) 

By introducing the integrating factor Yj^^ into (|1.45p . we obtain the following 
equation: 

The integration of Eq. (|1.45p yields 

= - ^^Y2 In \Y2\ + KY2 - (1.48) 
a 4a 

where K is an arbitrary constant determined from the initial conditions. 

Thus, expressions (|1.44p and (|1.48p determine a functional interrelation be- 
tween the price and the volume of the commodity, which is a permanent balance 
between the income of the producer and his expenses at any moment of time. 

Returning to the system (|1.38p . we should be reminded of the fact that the 
coefficients c?i, ^2 are the demand elasticity and its derivative with respect to 
the price, whereas Si, S2 are price demand elasticity and its derivative with 
respect to the volume of the commodity, i.e., 

di^7]D, d2^Tj'u, Si^rjs, S2 = v's- 

(It is understood that the elasticities and their derivatives are evaluated at 
the trivial state of equilibrium.) 

Then, conditions (|1.42p take a somewhat different form: 

Vd = b^Vs^ 

v'd = b'v's- (1-49) 

The first expression in (|1.49p implies a transition from a stable state of 
equilibrium to an unstable one with the formation (or annihilation) of a limit 
cycle, whereas the second expression in ()1.49p determines the boundary between 
the types of stability of self-oscillation regimes. 
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Chapter 2 



Periodic regimes in 
nonlinear models of the 
multiplier-accelerator 



2.1 A multiplier-accelerator model with finite 
duration of an investment lag 

In this Chapter, we consider an example of a certain kind of multipher-accelerator 
models that illustrates at a phenomenological level different aspects of the prob- 
lems of business cycles. As the starting point, we use the investment model 
of Goodwin [1] that possesses a nonlinear element built in the system of a 
multiplier-accelerator interaction. The dynamics of this model is characterized 
by the lag of two types: from the point of view of investment demand, there is 
a finite-duration lag of the action of the accelerator, whereas there is a contin- 
uously distributed lag on the part of supply. 

As the main variable, the model involves profit or an output Y = Y (t). An 
excess demand is formalized by the equation 

D = C + I + G, (2.1) 

where _D is a cumulative demand comprising a consumer demand C, an invest- 
ment demand /, and independent expenses G. All the above-mentioned terms 
are the actual costs. It is assumed that the consumption C is directly pro- 
portional to the profit and the lag is absent, i.e., C = cY, where the quantity 

< c < 1 characterizes a marginal propensity to consume. We also assume 
that the actual investments outlay is carried out with a certain fixed lag of 
the duration of 6 units of time after an investment decision F (t) is made, i.e., 

1 {t) = F [t — 0). Exactly here, the action of the accelerator manifests itself as 
a functional relation between the volume of investment decisions F (t) and an 
instantaneous velocity of the change in the profit (output volume) . In the 
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most general form, this relation can be represented as follows: 

fdY {ty 



F(t) = ip 



\ dt 



where is a certain nonlinear function possessing the property of saturation. 

In other words, for small changes in the profit, the classical linear accelerator 
comes into play, whereas for a further increase in the profit (output volume), 
the function F reaches its upper bound -Fmax determined by the resource and 
capacity limits of the structure of the production. 

When the volume of the production output substantially decreases, the quan- 
tity F tends to its lower bound i^min that depends, generally speaking, on the 
amortization quota of the fixed capital. For the sake of convenience of further 
mathematical transformations, within the framework of the model under con- 
sideration, we restrict ourselves to a Taylor expansion of the function Lp {x) up 
to the third order: 

ip {x) w -p^x -f (/?2 Y + (/^sy . 

where i^j, i = 1, 3, are corresponding derivatives of the function Lp (x). 

Given that we have already described the main components of the cumulative 
demand function, we obtain the following equation: 

D{t)^cY{t)+v{^^^+G. (2.2) 

Concerning the independent expenses G, it should be noted that they can 
be considered fixed, i.e., G — const. 

Consider next the situation on the part of supply. The main assumption 
is that value of the output volume Y lags behind the value of the cumulative 
demand D. The lag is considered to be continuously distributed, and it can be 
represented by a linear first-order differential equation: 

e'l^D-Y, (2.3) 

where e > is a constant of the time lag characterizing the dynamics of adjust- 
ment between the demand and the supply. 

A balance synthesis of the demand and the supply, by (|2.2|) and (|2.3p , yields 
a differential equation with a retarded argument: 

^^-(^)-<'-^' — »■ <-> 

Equation (|2.4p represents the general multiplier- accelerator model of Good- 
win with nonlinear interactions. 

With regard to the form of (|2.4p . it is nothing but a mixed differential- 
difference equation. 
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As is obvious, equation (|2.4p has a singular solution Y* = that represents 
an equilibrium level of profit (output volume) resulting from the action of the 
static multiplier. In what follows, it is reasonable to introduce the variable 
y (0 — ^1*-^'^ that represents a deviation of the profit from its equilibrium 
value. 

In terms of the new variable y (t) , equation (|2.4p takes the form 

ss'-^-Js'-y^)+s^yit)^0, (2.5) 



dt \ dt ^ 

where s = 1 — c is a marginal propensity to save. 

The given model, represented by Eq. (|2.5p . determines the dynamics of profit 
(output volume) in terms of deviations from the above-mentioned equilibrium 
level Y* ^ -r^. 

1 — c 

For the sake of a further analysis of dynamic properties of the considered 
process of changes in profit, it is necessary to carry out a sequence of mathemat- 
ical transformations and simplifications that will allow us to reduced the mixed 
first-order differential-difference equation to an ordinary differential equation of 
higher order [35j . 

Let us introduce new time t = t — 9 instead of t. Then, equation p.Sp can 
be represented as 

es^JlSZ±^ _ ^ (/JIP\ + ^r + e)= 0. (2.6) 



dr \ dr ^ 

Our next step consists in expanding the left-hand side of Eq. (j2.6p into a 
power series in 9, while retaining terms containing the first power of 9. We 
obtain: 

dyiT-) .<Py{T) of .dy(T)\ f dy(T)\ „ 



or 



In Eq. (|2.7p . we substitute the explicit form of the function ip as a. cubic 
polynomial. As a result of transformations, we obtain an ordinary second-order 
differential equation: 

d'^y ^ (£ + 9s~ip^\dy ^ sip^ fdyV s^Vz fdy\^ 



dr^ V J dr e9" 2e9 \dT J 6e9 \dT ^ 

Making the change of variables y = yi, ^ = 1/2, we represent (|2.8p in the 
form of a system of two differential equations: 

dyi 
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dy2 S Lp^-e-Os SLP2 2 , 3 lo n\ 

^ = --0^' + — re — + 2^^^ + fe^^- 

It is natural to begin the qualitative study of the dynamic system (j2.9p from 
an investigation into the states of equilibrium. As is obvious, equations (|2.9p 
possess only the trivial state of equilibrium yl = 0, t/2 = 0- To classify the type 
of this singular point, it is necessary to find characteristic numbers of the linear 
part of (|2.9p determined by the quadratic equation 

From the explicit form of (|2.10p . one can infer that the singular point of 
the system (|2.9p may be either a stable (unstable) node or a stable (unstable) 
focus. Of primary interest for us is the situation when a complex focus changes 
its stability, which may be accompanied by the formation of a limit cycle giving 
rise to a corresponding self-oscillation regime. In this case, we represent the 
solution of (|2.10p in the form 

Ai,2 = ^Ti^, (2.11) 

where = jg, = — 1, and /Lt = '^^ is a small quantity. In other words, 

for /i = 0, accordingly, the eigenvalues are purely imaginary: Ai^2 = ^iuj. 

By differentiation expression (|2.10p with respect to the parameter /x, we 
obtain: ~ ^ — \- This means that the eigenvalues Ai_2 cross the imaginary 
axis with a nonzero velocity. Therefore, we can argue that the conditions of 
Hopf's bifurcation theorem are fulfilled, and the system (j2.9p allows for the 
formation of a limit cycle from the complex focus. 

It would be in order here to draw attention to the reason for the occur- 
rence of instability in the multiplier-accelerator model. As it seems, exactly the 
accelerator "blows up" the damped oscillations induced by the multiplier and 
generates a structural self-sustained oscillation motion (self-oscillations). We 
have already seen that instability occurs in a given economic system when one 
of its parameters changes. It is most natural to consider as a variable parame- 
ter the coefficient of the linear accelerator (f^ whose critical value (^* = e + 9s 
changes the direction of damping in the system (|2.9p . Thus, the coefficient ip-y 
plays the role of a bifurcation parameter. The criterion of stability can now be 
written down in the form Lp^ < Lp\. A linear analysis shows that when tp^, while 
increasing, passes through the value ip\ , a loss of stability of the focus is caused 
by a pair of complex conjugate eigenvalues of the matrix of the coefficients of 
the linear part of the system (|2.9p . 

The loss of stability &i ip^ = ip^ occurs under the conditions of Hopf's theo- 
rem that states that, in addition to a stationary solution, there appear periodic 
solutions. 

However, Hopf's theorem itself does not provide information on whether 
these periodic solutions describe regimes that can be actually observed as steady 
ones. Periodic solutions may prove to be unstable and, accordingly, unobserv- 
able without the use of special procedures. Therefore, the next goal of our study 
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of emerging periodic solutions in the system (j2.9p is to derive explicit formulas 
describing their stability, amplitude, and period. 

To achieve the above-mentioned goal, we shall use the techniques presented 
in [33. For the sake of convenience of the application of the proposed methods, 
we shall retain original notation. 

In order to reduce the system (j2.9p to the normal Poincare form, we make 
the change of variables yi = xi, y2 — ~ujX2. As a result, for /i = 0, we obtain 

ii = -UJX2, 

2 , ^"^5^3 3 to^O\ 

X2=Loxi 7;—X2^ -T. — ^2- (2-12) 

2 D 

Let us represent the system of the two ordinary differential equations (|2.12p 
in the form of a complex differential equation with respect to the variable z — 
Xi + X2: 

z = iujz + g2o— + giizz + go2— +530 "7" +521^ 
2 2 D 2 

zz'^ z^ 

+ 512— + 503 (2.13) 

where z is the complex conjugate of z, and the inverse change of variables yields 
xi = Rez, X2 = Imz. 

For further evaluation, we need only the coefficients 520, 511, 502, and 521. 
Their explicit forms are 

520 = -.gii = .go2 = ^ , .g2i = • (2.14) 

Given expressions (|2.14p . we obtain the following: 
1) the value of the first Lyapunov quantity 

5^3 



16£26l 

2) the amplitude of small oscillations 



Reci (0) = (2.15) 



-5^3 



3) an approximate value of the period of oscillations 

27r A , sVi 



(2.16) 



T«— 1 + . (2.17) 

The periodic solution itself, up to the choice of the initial phase, in terms of 
the original variable, is written down in the form 
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It is important for us to know the sign of the coefficients ip^ that determines 
the sign of the first Lyapunov quantity in (|2.15[) . For ip^ > 0, Reci (0) > 0, 
and, accordingly, the Hmit cycle is unstable; that is, rigid excitation of self- 
oscillations takes place, accompanied by the phenomenon of hysteresis. How- 
ever, the assumption of the positivity of tp^ is unrealistic, because the condition 
of achieving the limit saturation Fmax (-Fmin) under an increase (decrease) in 
the velocity of changes of profit Y will not be satisfied. Therefore, we should 
set (fi^ < (Reci (0) < 0). In this case, a stable limit cycle is generated with a 
soft excitation regime of self-oscillations. 

While analyzing the explicit form of the approximate solution for Y (t) in 
expression (I2.18p . we should note the contribution of the coefficient (p2- A 
nonzero value of introduces certain asymmetry into the structure of the 
resulting oscillations. Obviously, they are nonharmonic even for small values of 
the amplitude. Besides, the coefficient induces an increase in the period of 
oscillations with a growth in their amplitude. 

An economic meaning of the asymmetry of the cycle consists in a difference 
between the duration of periods of expansion and that of periods of decline, 
which, on the whole, is a characteristic feature of nonlinear models of economic 
dynamics. 

Somewhat earlier, we considered in detail the influence of the accelerator 
parameter ip-^ on the degeneracy of the linear part of the system p.9p that 
directly induced a bifurcation of limit-cycle generation from an equilibrium state 
of the complex-focus type and the establishment of the regime of self-sustaining 
oscillations 12J. 

Summarizing, we want to emphasize that, in the present study of behavioral 
properties of the nonlinear multiplier-accelerator model of Goodwin, we have 
ascertained the mechanism of the occurrence of a cycle; we have determined 
the type of its stability, and we have given meaningful interpretation of the 
influence of all the parameters of the nonlinear accelerator on peculiarities of 
the self-oscillation regime. 

2.2 A model of the multiplier-accelerator with 
a continuously distributed lag 

In the previous Section, we have studied in detail an example of one nonlinear 
model of the multiplier-accelerator characterized by the occurrence of a stable 
limit cycle in the neighborhood of an equilibrium state. In that case, an essential 
element in the construction of the accelerator has been a finite-dimension lag 
in the functional relation between the measure of the volume of investment 
decisions and the instantaneous velocity of the change in profit (output volume) . 
As basic assumptions in the synthesis of the initial model, we have employed 
nonlinear dependence of the accelerator on the derivative of profit, as well as 
linear dependence of the consumption function on the value of profit. 

In the present consideration, the multiplier-accelerator model will be rep- 
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resented in a different form. First of all, we assume that the action of the 
accelerator is now expressed in terms of a continuously distributed lag [14] . 

According to R. Allen [1] , if the investment function / (t) represents the 
actually induced investment at the moment of time t, caused by changes in the 
output volume Y (t) , the lag is described by the differential equation 

where ip is a, nonlinear accelerator function, and the parameter f3 > character- 
izes the rate of changes in the investment function I (t). Concerning ip (^), it 
should be noted that its behavior for small changes in the profit is close to linear 
one, whereas with a further increase in ^ the accelerator function grows slower 
and may even become non-monotonic. We shall assume that the mechanism of 
the action of the accelerator is satisfactorily described by the cubic parabola 

V (u) = Piu + </?2 Y + <y53y > 

where (^j, i = 1,3, are corresponding derivatives of Lp [u). 

The next step consists in the introduction of a lag into the model of the 
multiplier. In an analysis that follows, we shall also employ a continuous repre- 
sentation for the description of the action of the multiplier effect by means of a 
corresponding differential equation. 

Let us assume that on the part of cumulative demand, a lag is absent. The 
planned consumption is given by C = C(y), and independent expenses are 
determined by the quantity G. 

Then, the cumulative demand can be represented in the form of the equation 

D^C{Y)+I + G. (2.20) 

Here, we give up the hypothesis that the consumption function is linear and 
assume that C (Y) is a substantially nonlinear function of the output volume 
Y. 

In the most general case, we restrict ourselves to representing C (Y) in the 
form of a third-order polynomial, i.e., 

C{Y)^CiY + C2— + C:i — . 

The coefficients Ci, i = 1,3 have the same meaning as in the case of the 
accelerator function. Moreover, Ci > 0, whereas Ci and C2 may have opposite 
signs. 

The parameter of independent expenses G is considered to be constant. 

Consider the situation on the part of supply. Here, a response of the out- 
put volume Y to the cumulative demand is considered to be non-instantaneous, 
inertial, i.e., there exists a continuously distributed lag in the form of a corre- 
sponding differential equation: 

^^aiD-Y), (2.21) 
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where ^ is a time constant of the muhipher 

Equations (|2.19p - (|2.2ip completely determine the system of two nonlinear 
differential equations that describe an interaction between the multiplier and 
the accelerator. This system is represented as follows: 

dY 

— ^aiC{Y)-Y + I + G), 



dt \ \ dt ^ 

Given that <^ (0) = 0, the system (I2.22[) has singular solutions determined 
by the system of algebraic equations 

C{Y)-Y + G = Q, 

1 = 0. (2.23) 

The first equation of (|2.23[) , by the explicit form of C (Y) , may have up to 
three roots that represent the coordinates of an equilibrium state of the system 

We represent the system of two differential equations (|2.22p in the form of 
a single second-order differential equation: 



fY_ 



(2.24) 
where C" (F) = 

With the help of (|2.24|) . we have succeeded in eliminating dependence on 
the variable I{t). In what follows, we shall operate a vector field composed 
of the variables Y {t) and '^^Jf '' ■ Introducing the coordinates yi = Y (t) and 
2/2 = ^^^^t^i it is reasonable to transform Eq. (|2.24p into a system of two first- 
order differential equations: 

yi = y2, 

y^ = a(^(^C'{yi)~l-^^y2 + P{v{y2) + C{yi)-y, + G)^ . (2.25) 

As is obvious, the system (|2.25p possesses the same states of equilibrium as 
the system (|2.24p does. 

Let the coordinates of a single point be given as (2/1,2/2 — 0), where 2/1 is 
the solution of the equation C (yi) ~ yi + G = 0. In order to analyze behav- 
ioral properties of the system (|2.25p in the neighborhood of the given state of 
equilibrium, we introduce new variables ui = yi — y* and U2 =2/2- 

Given the explicit form of the nonlinear consumption function C iY) , as a 
result of some transformations, we arrive at the system of differential equations 



ill = U2, 
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V? 

U2 = -al3Sui + - (3 - aS) U2 + a/3 {C^yl + C2) 



1 



+a {C-iy\ + C2) U1U2 + a(3(p2^ + a/JCs^ 

2 

+ aC3^ + aMf , (2.26) 

where 5 = 1 - Ci - C22/1* - Cs^- 

Let us consider a particular case of the system (|2.26p under the conditions 
(P2 = 0, C32/* + C2 = 0. Under the above-mentioned restrictions, the quadratic 
terms in (|2.26p vanish; that is, there is symmetry with respect to the change 
ui <-> — ui, U2 ^ — U2. Furthermore, under the assumption that 

C2(C| + 3C3(l-Ci)) 

the equation for singular points j/i is factorizcd as follows: 

(C32/1 + C2) + —2/1 j = 0. (2.27) 

Accordingly, yl = is the coordinate of the state of equilibrium in whose 
neighborhood the behavior of the system (|2.26p is being studied. To ensure 



the positivity of yl, we assume that C2 < 0, C3 > 0, and the parameter 5' = 
l-Ci + ||>0. 

After the above-mentioned simplifications, the system (|2.26[) takes the form 



Ul = U2, 



U2 = —afSSui + (aPipi — P — aS) U2 + afiC^ — 

6 

+ aC3^-)-a/3^3|. (2.28) 

We shall be concerned with qualitative properties of the system (j2.28p in the 
neighborhood of the trivial state of equilibrium u\ = 0, U2 = 0. To determine 
the type of equilibrium, it is necessary to ascertain spectral properties of the 
linear part of (|2.28p with the characteristic equation 

- /iiA -I- a/JS" = 0, (2.29) 

where Hi = a(3(pi — aS — (3. 

The explicit form of the quadratic equation ()2.29p is analogous to that of 
the characteristic polynomial in [12 . Therefore, when analyzing the situation 
with the occurrence of a periodic regime in the system (I2.29P resulting from the 
change of stability of the singular point of the type of a complex focus 
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we can arrive at the conclusion that the reason for this effect is transition of the 
linear parameter of the accelerator through a certain critical value (^f — "^^^ . 

Note that the derivative of the eigenvalue A with respect to the parameter 
fii is nonzero: 

dX _ 1 _/ Q 
dfi^ 2 

In this case, we may argue that the conditions of Hopf 's theorem are satisfied, 
and a limit cycle around the trivial state of equilibrium is generated in the 
system (I2.28|) from a complex focus. 

The fact that a self-oscillation regime is present is rather remarkable in itself; 
however, it does not provide much information. Hopf 's bifurcation theorem does 
not give any answer to the question about the uniqueness of the limit cycle and 
the character of its stability. 

To resolve the posed problems, we reduce ()2.28|) to the normal Poincare form 
for fi^ — 0, using the change of variables ui — xi, U2 = —ujxi: 

X2 = UJX2 — + aCa— 1- aptp^uj —. (2.30) 

Lu b 2 6 

The system (|2.30l) can be reduced to the complex differential equation 

Z3 Z^Z ZZ^ Z3 , , 

Z = iujZ + .930-^ + .921—^ ^^^^~2 ^903-^, (2.31) 

where Z = xi + 1x2, Z = xi — 1x2] 

a (3C3 - /3C.V3 - ^) ^ _ « (C3 + /3C.V3 - ^) ^ 



330 — z ' 321 



912 — ^ ; .903 — ^ 

Making use of the explicit expressions for the coefficients gij , it is not difficult 
to determine the first Lyapunov quantity: 

/,(0) = Ref = (2.32) 

As before [5], we assume that, as a result of the effect of investment sat- 
uration, < 0, whereas the coefficient C3 > 0. Therefore, a sign change 
in expression (|2.32p is possible, which is a manifestation of different types of 
stability of the limit cycles. From (|2.32p . the stability of the limit cycle for 
C3 -I- f3uj^(f^ < {li (0) < 0) follows directly, whereas for the opposite sign of 
the inequalities an unstable self-oscillation regime with a catastrophic loss of 
stability takes place. 
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The case when the first Lyapunov quantity is small and alternates the sign, 

i.e., 

h = /i2, (2.33) 

is of much greater interest. 

As is well-known [6], the behavior of dynamic systems in the vicinity of 
the parameter values satisfying the equality C3 + (iuj'^'p^ = 0, such that the 
first Lyapunov quantity h vanishes, substantially depend on the sign of the 
second Lyapunov quantity I2. Depending on the first and the second Lyapunov 
quantities, as well as on the sign of the real part of the roots of the characteristic 
equation /i^, in small neighborhood of the state of equilibrium on the phase 
plane, one or two limit cycles may exist with all possible combinations of stability 
and instability: namely, one stable or unstable limit cycle, or two limit cycles 
(a stable one inside an unstable one or vice versa). 

The second Lyapunov quantity is determined by the expression 

^2 = -7^Im.g3o5i2- (2.34) 
After substitution in (|2.34p of the actual values of the parameters, we obtain: 

or, by the vahdity of C3 — Poj^(p3 — 0, 



3^2 



(2.35) 



As is obvious, the quantity I2 does not vanish for any values of the parameters 
and is strictly positive, i.e., I2 > 0. 

If we make a conversion from the complex- valued variables to polar coordi- 
nates, we get two independent equations for the amplitude and the phase of the 
cycles: 

P = P (Mi + M2P^ + hp'^) , 

^ = uj. (2.36) 

The states of equilibrium for the first equation of (|2.36p satisfy the bi- 
quadratic equation 

Ml + /^2P^ + ^2P^ = 0. (2.37) 

Equation (j2.37p may have either none or one, or two positive solutions (cy- 
cles). 

In Fig. 2.1, we present the corresponding bifurcation diagram. The line 
H = {(M11M2) ■ Ml — 0} relates to the usual Hopf's bifurcation. The state of 
equilibrium is stable for < 0, and it is unstable for fii > 0. If we move along 
the line fx^ — to the points where < 0, the complex second-order focus 
on the phase plane will generate an unstable (coarse) limit cycle, whereas the 
focus itself becomes non-coarse and stable. Should we enter, while crossing H~ , 
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1^2 

H 




1^1 



Figure 2.1: The bifurcation diagram. 

region 2, the stable complex focus generates a stable limit cycle. In region 2, 
both the cycle, the stable one and the unstable one, coexist simultaneously; the 
merge and disappear on the line T = {(/^i,/J.2) : A*2 = 4/2/^^, ^2 < O}- 

The line T characterizes the bifurcation of the double cycle. Further in 
region 3, limit cycles are absent. 

In Fig. 2.2, the region of the coexistence of the two limit cycles is shown. 

In conclusion of this investigation, we would like to point out the following: 
the presence of two limit cycles in the initial dynamic system of the multiplier- 
accelerator is stipulated not only by the nonlinearity of the accelerator function, 
but also by substantial nonlinearity of the consumption function, because ex- 
actly a relation between nonlinear coefficients of these functions ensures the 
bifurcation of the double cycle. 

2.3 Cyclic regimes in a nonlinear model of the 
multiplier-accelerator with two degrees of 
freedom 

Consider a model of the multiplier-accelerator with spatial inhomogeneity. Such 
a model reflects peculiarities of interregional trade in the presence of an import- 
export multiplier, which agrees with the already studied multiplier of local ex- 
penses, as well as with the general concept of Keynes and Samuelson [30l [31] . 
Let the quantity of an imported commodity depend on the local profit which 
is now a function of time t and of a generalized spatial coordinate r. Assuming, 
as the first approximation, that the action is local, we suppose that the commod- 
ity is imported from the nearest neighborhood of a considered point, whereas 
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Figure 2.2: The stable and unstable limit cycles. 



the exported product is produced under the influence of the same propensity to 
import in the neighborhood of this point. Then, the net trade surplus is deter- 
mined by the product of a constant propensity to import and a profit margin 
in export-import operations. In other words, the above-said can be represented 
in the form of the following expression: 

X — M — m— ^, 

where X is an export volume; 
M is an import volume; 
m is a constant propensity to import. 

The spatial model of the multiplier-accelerator is represented in the form of 
a single second-order differential equation with a nonlinear investment function: 

q2y q2y dY V /dY\^ 

TnT^Sy-.^^ = W-l-S)^--(.^) , (2.38) 

where 5 is a marginal propensity to save; 
1/ is the coefficient of the accelerator. 

In our consideration, we assume that all the parameters of the model, i.e., 
5, u, and m, are constant positive quantities that do not depend either on time 
or on the spatial coordinate. 

The model described by Eq. (|2.38p is a rather complicated mathematical 
object that exhibits a variety of forms of spatial-temporal organization. There- 
fore, in what follows we shall focus on dynamic processes, having preliminarily 
subdivided the space into two parts interrelated by the regional trade. This will 
allow us to study such phenomena as frequency matching and quasi-periodic 
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motion. In other words, we may encounter a new form of the attractor, namely, 
an invariant torus. 

The partial differential equation (|2.38p is represented in the form of two 
coupled ordinary differential equations of the second order: 

Fi + {Si + mi) Fi - m2Y2 ^ (ivi - I - Si) Yi - '^Y^^ 

% + (^2 + ms) Y2 ~ miYi = {v2-l- S^) % - y^a'- (2-39) 

Here, we have assumed that the parameters of the accelerator, the rates of 
accumulation and of import are different for each region of the subdivision. 

Previously, for the model of the multiplier-accelerator in one region, we have 
found periodic regimes with the emergence of corresponding self-oscillations, 
and we have studied the character of their stability. As it seems, for the model 
with two regions, cyclic motion is also possible. Moreover, quasi-periodic motion 
with two matched frequencies is possible as well. 

Let us represent the system (j2.39p in the traditional form of a system of 
first-order differential equations. This new system is four-dimensional. 

Let x\ = Fi, X2 = Yi, X3 = I2, and X4 = Y2- As a result, the system (|2.39p 
takes the form 



Xi = X2, 



2 



±2 = - {Si + mi)xi + (j^i - 1 - Si) X2 + m2X3 - ^^i 3 . 
X3 = X4, 

X4 = mixi - (5*2 + m2) 2:3 -I- (!^2 - 1 - 5*2) X4 - V2-^- (2.40) 

Obviously, the system (|2.40p has trivial equilibrium x* = 0, i = 1,4. The 
matrix of the linear part of (|2.40p . corresponding to this singular point, has the 
form 

{010 
-(5i+mi) z^i-l-S*! m2 ^ ■ , . 

^ }■ i2.41) 

mi -(5*2+7712) 

The matrix (|2.4ip has the characteristic polynomial 

(A' -{vi-l-Si)\ + Si+ mi) (A^ 

- (1^2 - 1 - ^2) A + 52 + m2) - mim2 = 0. (2.42) 

Equation (|2.42p is a fourth-order equation; hence, it has four roots. Of 
primary interest for us is the situation when (|2.42p has two pairs of complex 
conjugate roots with small parameters in their real parts, i.e., 

(y.\ 0.2 
Ai,2 = — ±iwi, \z^4: = — ±iL02- (2.43) 
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Under the assumption that = I + Si + ai, i = 1,2, there exist, for 
ai = 0, critical values of the parameters of the accelerator i>'^ = 1 + Si that are 
responsible for a possible formation of limit cycles. 

For ai — a2 = 0, equation ()2.42|) reduces to to the biquadratic equation 

A"^ + (Si + mi + 5*2 + m2) + {Si + mi) {S2 + 7712) - mim2 = 0, 

which, for = ^t^^, yields an equation for the frequencies: 

uj'^ - {Si + mi + 5*2 + m2) + {Si + mi) (^2 + m2) - mim2 = 0. (2.44) 

As the free term in ()2.44p is positive, equation (I2.44|) has two positive roots 
that determine the squared frequencies: 



D — {Si + mi + S2 + m2)^ — 4 (5i + mi) (^2 + m2) — mim2 > 0. 

For definiteness, we assume that uji > lOi- If we compare the values of the 
so-called eigenfrequencies (D^ — \J Si + mi, i — 1,2, it is not difficult to prove 
that uJi > uJi, whereas L1J2 < uj2- It means that the common frequency of the 
matched oscillations is either higher than the maximum natural frequency or 
lower than the minimum natural frequency. On the other hand, in the case 
of free oscillations, the matched system will never be able to oscillate with an 
intermediate, compared to natural frequencies, frequency. From an economic 
point of view, this fact means that a connection between two regions by means 
of trade relations either speeds up or slows down a cycle of business activity in 
both the regions. 

For further consideration of the properties of the four-dimensional flux that 
possesses a state of equilibrium with two pairs of purely imaginary eigenval- 
ues, it is necessary to construct the normal form for the system of ordinary 
differential equations ()2.40p . This can be done by means of a sequence of linear 
transformations of the initial variables as follows: 




(5i+mi + 5'2+m2)±Vi^ 
2 



xi = m2 {yi + ys) , 



X2 = -m2 {uJiy2 + ^22/4) , 





The transformation (|2.45p converts the system of four ordinary differential 
equations in real variables into a system of two complex differential equations 
that takes the following form in terms of the polar coordinates Zj = p^e^'^^: 



Pi = CKiPi + aiipl + ai2Pipl, 



P2 = "2^2 + 0'2ipIp2 + a22pl 



36 CHAPTER 2. PERIODIC REGIMES IN NONLINEAR MODELS OF THE MULTIPLIER-ACCELE, 



Lp^ = UJl 



LP2—UJ2 + 



o{\p\'). 



(2.46) 



Here, \p\^ — Pi ^ P21 "^ii are small sign-alternating parameters; the 
coefficients Zj ~ p^e^'^j , j — 1,2, are functions of the initial parameters of the 
system. 

We can learn a lot about the dynamics of the system (|2.45p from the consid- 
eration of a plane vector field obtained by discarding the angular coordinates, 
following the methods proposed in [42] . 

In order to reduce the number of parameters, we scale the variables pi and P2- 
Setting pi = pi ^y\all\ and P2 = P2 v|a22| , dropping for notation convenience in 
what follows the bar over p^ and P2, and, if necessary, scaling the time variable, 
we obtain: 

Pi = Pi ("1 + pI + bpl) , 



P2 = P2 ("2 



cpl 



dpi) 



(2.47) 



±1, 



ai2 
1 022 1 



Q21 
|aii| 



The system (I2.47|) is characterized by twelve topologically different situa- 
tions, presented in the following table: 



Case 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


d 


+ 


+ 


+ 


+ 


+ 


+ 














b 


+ 


+ 


+ 








+ 


+ 


+ 








c 


+ 


+ 




+ 






+ 






+ 


+ 




d — be 


+ 




+ 


+ 


+ 






+ 




+ 







This classification is based on an analysis of secondary "pitchfork" bifurca- 
tions from nontrivial states of equilibrium of the plane vector field. Note that 
the singular point {pi, P2) — (0,0) is always a state of equilibrium; besides, up 
to three states of equilibrium may exist in the positive quadrant: 



{Pi > P2) = (V-ai, 0) , for ai < 0; 
(Pi)P2)== 0, W— ^ , forQ!2d<0; 



(2.48) 



iPl,P2) 



I da2 — bai / cai 



a2 



Q 



for 



da2 — bcti cai — a2 



Q 



Q 



<o. 



where Q — d — be, d ~ ±1. 

In general, the behavior of the system remains comparatively simple until 
the occurrence of secondary Hopf bifurcations from the fixed point {pi , P2) — 



Q 



^"^Q ). In order to detect such bifurcations, we linearize the 
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dynamic equations in a small neighborhood of this singular point. As a result, 
we obtain the following matrix: 



whose trace is 



trr= ^(aid(c-l) + a2 {b - d)) . 



and the determinant is 

4 

detT — — {{ba2 ~ da\) (cai — 02)) . 

Taking into account the conditions of the existence of the singular point 
(|2.48p . we infer that the secondary Hopf bifurcation may occur only on the 
straight line 

"2 = (2-49) 
— a 

and, at that, Q > 0. 

From this fact, it follows immediately that the secondary bifurcation does not 
occur in cases 2, 6, 7, 9, 11, 12. It is also possible to show that this bifurcation 
does not occur in cases 1, 3, 4, 5, because, for its realization, it is important 
that the angular coefhcient of the straight line (|2.49p should lie in between the 
angular coefRcients of the straight lines that correspond to the "pitchforks", i.e., 

d 

a2 = cai, a2 — t<^1: (2.50) 




which is equivalent to 



d(l-c) d 



in a corresponding sector of the plane (ai, (22 )■ 

As can be shown by means of simple evaluation, in each case, this require- 
ment does not ensure the condition Q > 0. 

Consider case 8, in which a Hopf bifurcation may occur. Some bifurcation 
sets and phase portraits for this case are represented in Fig. 2.3. 

On the Hopf-bifurcation line (|2.49[) . the system 

Pi = Pi ("1 +pI + bpl) , 

is integrable, whereas the function 

^(Pi,P2)-PiP2 ("i+P?+7P2), (2.52) 
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with 9 — P = ^'^d 7 = is constant along the solutions. In 

case 8, we have: & > > c, Q = — 1 — 6c > 0, and ai = —a < 0; therefore, level 
lines of this function have the form shown in Fig. 2.4. 

As the system (I2.5ip is integrable, the secondary Hopf bifurcation is de- 
generate. Therefore, to study the topology of this bifurcation in full detail, it 
is necessary to consider in (|2.49|) terms containing higher powers of the initial 
variables. 

Unfortunately, we do not possess necessary information, because nonlinearity 
in the accelerator model is restricted to cubic terms, whereas an exhaustive 
analysis requires terms of the fifth order. Nonetheless, it is possible to discuss 
some conclusions that can be drawn on the basis of the obtained results about 
the total four-dimensional flux of the system (|2.46p . It is necessary to bear in 
mind that here we have two rotations ifi — lui and (/?2 = Moreover, the 
ratio of the frequencies uii : uj2 that should be restored for final conclusions. As 
is easy to see, the initial four-dimensional system (|2.46p possesses four types of 
attractors corresponding to the fixed points of the plane system (|2.47l) : 

1) the trivial fixed point {pi =0,p2 = 0); 

2) a periodic orbit with the period ~ ^ (Pi = Pi, P2 = 0); 

3) a periodic orbit with the period ~ ^ (Pi = 0, P2 = P2); 

4) an invariant two-dimensional torus with the periods 

« — , « — (Pi = o,P2 = P2) ; 

5) an invariant three-dimensional torus with the periods 

uji 0J2 \a J 

The last (long) period on the tree-dimensional torus is associated with a 
secondary Hopf bifurcation in the plane system (|2.47p . Here, we can foresee 
the presence of subtle resonance effects and phenomena that remain outside the 
scope of our consideration. Thus, we expect to find a narrow "wedge" around 
the line secondary Hopf bifurcation where chaotic dynamics, including transver- 
sal homoclinic orbits and " horseshoes" , takes place. Consequently, we can argue 
that the multiplier-accelerator model of business cycles with a nonlinear invest- 
ment function, extended to the case of interregional trade, and whose motion is 
induced by a linear multiplier of import, may initiate random motion. 

This fact implies that processes of economic forecasting are problematic. As 
it seems, short-term forecasts are the most efficient ones, because exponential 
divergence of close trajectories does not occur over small periods of time. 
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Chapter 3 



Self-organization in 
Keynesian models 

As a science, macroeconomics was born owing to J. M. Keynes and, in the first 
place, to his outstanding work "The General Theory of Employment, Interest 
and Money" where, for the first time, the problems of macroeconomics appeared 
as the main subject of the research. Keynes's viewpoints, revolutionary at that 
time, had substantially formed as a result of an analysis of the reasons for the 
occurrence of the Great Depression. The global crisis of 1929-1933 prompted 
Keynes, as well as many other economists after him, to study seriously the 
economy as an integral system |27j . 

In this Chapter, we shall pay attention to some problems of economic dy- 
namics, the models of which are based on the postulates of Keynes's theory. As 
such, the following models will be successively considered: 

1) the model of the growth of the gross domestic product; 

2) the LS-LM Keynes model; 

3) the Kaldor model. 

3.1 The dynamics of GDP growth 

In the construction of long-term programs of social-economic development and 
in macroeconomic modeling, the most important, in terms of criteria, factor is 
the volume of GDP per capita. Following the post-Keynesian tradition, con- 
sider, as the basis of the model of GDP growth, the main macroeconomic identity 
for the volume of the total income: 

Y = C + I^G + Ne, (3.1) 

where Y is the GDP volume in value terms. 

Let us analyze the terms of expression (j3.ip from the point of view of their 
economic meaning. 
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1) C = C (y) is the consumption volume that depends on GDP. We shall 
assume that C (y) is a nonlinear function that can be represented by quadratic 
dependence: 

C (y) ^ Co + Ciy + C2y^ 

where Cq is autonomous consumption; Ci is the limit propensity to consume 
for small values of y; C2 characterizes extremal properties of the consumption 
curve. If C2 < 0, then C (y) has a maximum, which agrees with Keynes's 
postulate that the propensity to consume declines with an increase in income. 
For C2 < 0, the function C [y) has a minimum, according to the Modigliani- 
Duesenberry hypothesis of relative income [S^ . 

2) I = I {Y, ^) is the investment function. We assume that investment 
react only to the rate of changes in GDP, i.e., it carries out the function of a 
simple accelerator: 

dY 

^ = "^' 

where v is the marginal capital coefficient. 

3) G is the volume of government expenditure. In a simplified version, we 
assume that the quantity G is independent of GDP and is constant in time. 

4) Ne is the volume of pure export characterizing the factor of external 
economic activities. 

By analogy with G, we also assume that Ne = const. 

Taking into account all the above assumptions, we can represent the model 
(j3.ip in the form of an ordinary differential equation: 

dY 

v—^^R+{\^C^)Y^C2Y\ (3.2) 
dt 

where R — Cq + G + Ne = const. 

Singular solutions to Eq. (|3.2p that correspond to states of static equilibrium 
can be found from the condition 

or 

C2Y^ - {\-Ci)Y + R = Q. (3.3) 
The quadratic equation (|3.3p has the following representation for the roots: 

Y^*-2= ^~^2C2^ ' ^ = -4C2i^+(l-C•l)^ (3.4) 

If Ci < 1, which is usually assumed, and A > 0, then Y^2 ^-re different 
positive numbers. In other words, equation (|3.2p has two states of equilibrium, 
and the coefhcients C2 and R have the same sign at that. In the case A — 0, 
there exists a double state of equilibrium: Y^2 — "^202 ' 

For A < 0, equation (|3.3p has no real solutions, i.e., equation (|3.2p does not 
possess states of equilibrium. If the coefficients C2 and R have opposite signs. 
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one of the roots is a positive number, whereas the other one is negative. As 
the negative solution has no economic meaning in this case, one speaks about a 
single state of equilibrium. 

For the convenience of a further analysis of the properties of the differential 
equation (|3.2p , we introduce new variables X = Y — ' ^ ~ v' ^ result 
of the change of variables, equation (|3.2p takes the form 

Expression (j3.5p is a typical differential equation with separable variables. 
The integral curve that passes through the point t = 0, X (0) = Xq for C2 > 
is given by the following equations: 

^i^)= ,^r\ ^ if^ = 0; (3-6) 
1 + C2A0T 

4X0C2 + A tanh (VAt) 

X(t) = ^ ^, ifA>0; (3.7) 

4C2\/I + 4X0 C| tanh (VAt) 

4XoC2\/^ + Atan(\/^T) , ^ 

X(t) = %dl \ ^ \ , if^<0. 3.8 

It is not difficult to notice that, for different values of the quantity A, the 
solutions to the differential equation (j3.5p differ substantially with regard to 
their properties, i.e., even for small variations of A in the neighborhood of zero, 
a qualitative change in the scenario of the evolution of X (r) takes place. There- 
fore, the development of the situation can be diagnosed by means of qualitative 
methods, without resort to complicated and expensive calculations. In qual- 
itative forecasting, special attention should be paid to those factors that can 
change the dynamics of GDP growth either in negative or positive direction. 

In other words, if the process under consideration is in a zone of stable de- 
velopment, the rest of qualitative information (such as the actual trajectory of 
development) becomes less important. As a rule, in macroeconomic modeling, 
mistakes and errors of approximation in the parameters of the model are possi- 
ble; accounting for perturbations of exogenous character is also rather difficult. 
In this regard, it is reasonable to draw conclusions not about a single trajec- 
tory of development, but rather about a region of space of possible trajectories. 
The evaluation of such regions is the subject of qualitative theory of differential 
equations. 

Let us normalize the variable X (r) in such a way that will allow us to reduce 
the number of parameters in the differential equation (|3.5p to a single one. 
For X (t) ~ — ^^Tp- , equation (|3.5p takes the form 



^ = fiU,A), fiU,A) = U'-^. (3.9) 
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u 



A>0 



A = 



A<0 



Figure 3.1: The bifurcation diagram for the case of "double equihbrium" . 

Let U* = U* (A) be a state of equihbrium of Eq. (|3.9p for certain for certain 
fixed values of the parameter A, and let A (A) — f ([/*, A). 

For A < 0, the state of equilibrium U* is stable, whereas it is unstable for 
A > 0. For small variations of the parameter A, the behavior of the trajectory 
(|3.9P in the neighborhood of the state of equilibrium with A 7^ does not change 
qualitatively. 

Thus, the inequality A 7^ is the condition of non-degeneracy that singles 
out a coarse case. 

In the neighborhood of the state of equilibrium, the coarse system is modeled 
by the linearized equation p.9p : 



A different situation occurs when, for certain values of the parameter A, the 
eigenvalue A vanishes in the vicinity of the state of equilibrium: 



and the condition of non-degeneracy, f"uu 7^ 0, is fulfilled. In this case, U* is 
a double root of the equation / {U*) ~ 0. 

A model equation for this bifurcation depends on a single parameter and 
has the form (|3.9p . Then, for A > 0, the system (|3.9I) possesses two states 
of equilibrium: namely, a stable one and an unstable one. (The latter is the 
boundary of the attraction region of the stable state of equilibrium.) For A = 0, 
they merge into a "semi-stable" state of equilibrium and a non-coarse system 
appears. For A < 0, the states of equilibrium disappear: see Fig. 3.1. 

Let us monitor the stability of the state of equilibrium. Thus, when the 
parameter approaches the bifurcation value ^ = 0, the attraction region shrinks 
on one side, and, after the disappearance of equilibrium, all the solutions leave 
the considered phase space. In economic applications, such a phenomenon is 
called " a break-down of equilibrium" . 




(3.10) 



A(A) = /'(C/*,A)=0, 
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Note that, for the bifurcation value of the parameter A = 0, the projective 
mapping f{U,A) = onto the parameter space has a "fold" -type singularity. 
In other words, there occurs a loss of stability, i.e., a catastrophe of the "fold" 
type. 

An important quantitative parameter of the zone of stable development, in 
the model of economic development, is represented by the parameter A, which 
is a measure of the distance between the states of equilibrium Yi and Y2 (or Xi 
and X2): 

A = Ci (Yi - Y^f = Cl (Xi - X^f . 

In the case, when A is sufficiently far from zero, one of the states of equi- 
librium is stable, whereas the second one is unstable. Such a situation is rather 
typical of, for example, the study of logistic economic growth. When A tends to 
zero, the two states of equilibrium merge into a single one, and one can observe 
a state of dynamic uncertainty in the sense of the stability of economic growth 
trajectory, with a possible "ejection" of the corresponding trajectory into the 
region of depressive dynamics. Returning to the basic macroeconomic equation 
(|3.2p . we can formulate the condition of stability: the derivative of the con- 
sumption function with respect to the GDP volume at the point of equilibrium 
should be larger that unity: 

C'y{Y*)>1. (3.11) 

In this case, according to the Keynes absolute-income hypothesis, the con- 
sumption function is linear: 

C{Y) = Co + CiY. 

Then, from ()3.11|) . it follows that exponentially unstable growth is observed, 
because earlier we have assumed that the condition Ci < 1 is fulfilled. 

In other words, the model of economic growth with a linear consumption 
function can be stable only under the condition that a marginal propensity to 
consume is larger than unity. As is obvious, such a condition cannot be fulfilled 
in practice for an arbitrary GDP volume. This, in turn, induces a discussion 
about differences in the behavior of the consumption function in short-term and 
long-term periods. 

The consumption function we have mentioned before and which rests on the 
relative-income hypothesis is devoid of such shortcomings. 

Thus, the problem of studying the stability of economic growth, described 
by Eqs. p.ip and (|3.2p . under the above-mentioned substantial simplifying 
assumption, has been reduced to structural properties of the consumption func- 
tion. Quite naturally, there always exist certain parameters that, for some 
reason, do not fit into our formal analysis. This fact hampers the application 
of structural methods of macroeconomic forecasting, such as, e.g., simulative 
modeling. For this reason, in the present work, we employ the methodology of 
ascertaining, analyzing and forecasting economic processes and phenomena of 
self-organization by means of qualitative methods of the examination of dynamic 
models. By ascertaining substantially nonlinear processes of self-organization. 
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we are able to single out the controlled factors as well as those factors that are 
beyond out influence [34) . 

As a result of our investigation, we have clearly demonstrated that the 
macroeconomic model of growth may have several states of static equilibrium. 
Accordingly, in dynamics, there exist evolutionary trajectories that correspond 
to the states of equilibrium of the system. In the case when the structure of 
the economy is coarse, small perturbations of the environment are unable to 
"dislodge" it from its trajectory of growth. In the situation of a non-coarse 
system, even small jumps of the parameters (or a " shock" ) may cause a sudden 
transition to a different state of equilibrium, and the trajectories of development 
undergo qualitative changes. On the other hand, a merger of states of equilib- 
rium is accompanied by a bifurcation with a catastrophic loss of stability. In 
this work, for a quadratic consumption function, we have studied effects in the 
vicinity of a double state of equilibrium, with the occurrence of a "fold" -type 
catastrophe. Taking into account higher-order nonlinearity makes possible to 
observe a bifurcation of threefold degeneracy of the state of equilibrium, with an 
"assembly" -type catastrophe, of fourfold degeneracy of the state of equilibrium, 
with an "swallow-tail" -type catastrophe, etc. 

A qualitative study of properties of the macroeconomic model should precede 
the elaboration of efficient methods of the management of the economic system 
on the basis of self-organization principles. The use of self-organization allows 
one to optimize shortages of the functioning of the economic system and to 
prevent the control parameters from entering a zone of undesirable bifurcations 
and catastrophes [22]. 



3.2 The LS-LM Keynes model 

Consider the dynamic economic system proposed by J. M. Keynes in his work 
"The General Theory of Employment, Interest and Money" [20]. This model 
represents the conditions of mutual balance in the goods and money market and 
is called, accordingly, the LS-LM model. In its simplest form, a business cycle 
is described by the system of ordinary differential equations 

dY 

TY^ = HY,R)-S{Y,R), 

TR^=L(Y,R)~M. (3.12) 

According to W. Zang [19] , all the parameters and variables here are positive, 
and their meaning is the following: Y = Y (t) is the volume of the national 
income; R — R{t) is the interest rate; I = I (Y, R) is the investment demand 
function that increases with respect to the volume of the national income, i.e., 

= Iy > 0, and decreases with respect to the interest rate, i.e., ^ = —Ir < 
0; S = S (Y, R) is the savings function that increases with respect to both the 
variables, i.e., Jp- = Sy > 0, ^ = Sn > 0; L — L{Y, R) is the total demand 
for money that increases with respect to the income, i.e., ^ = Ly > 0, and 
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decreases with respect to the interest rate, i.e., ^ — —Lr < 0; M is the 
constant supply of money; Ty and tr are corresponding time constants. 

The system p. 121) illustrates the action of a simple mechanism: an excess 
demand for investment, compared to the volume of savings, leads to an increase 
in the national income, and vice versa; in the case when the total demand for 
money resources is higher than the available supply, the interest rate increases. 

Concerning the investment demand function, it should be noted that the 
values of / {Y, R) stand in a direct relationship to the volume of the national 
income and in an inverse relationship to the interest rate. On the other hand, 
this means that an increase in the income or in the interest rate will stimulate 
the population to enhance the savings; and under the condition of growth in 
the output (income) volume, the demand for money L (Y, R) increases. 

For the system p.l2p . we assume the existence of, as a minimum, one positive 
singular solution Yq, Rq, which represents the state of static equilibrium of the 
LS-LM model. For algebraic evaluation of the combination of the values of the 
national income and of the interest rate, it is necessary to solve the following 
system of equations: 

I{Yo,Ro) = S{Yo,Ro), 

L(Yo,Ro) = M. (3.13) 

It is sufficient to confine the analysis of dynamic properties of the system 
(|3.12p to a local domain of the two-dimensional space of the initial variables 
Y (t) and R (t) near the state of equilibrium Yq, Rq. 

To this end, we introduce new variables Y (t) = Y (t) — Yq and R (t) = 
R{t) — Rq, with the meaning of deviations from the equilibrium values of the 
income and of the interest rate, and expand the right-hand sides of the system 
(|3.12p in a Taylor series at the point of equilibrium, retaining the first and 
the second powers of the corresponding variables. For convenience, we drop 
the bar over the variables and, without loss of generality, set tr = ty = 1, 
F{Y,R) = L{Y,R)-S{Y,R). 

As a result, the system (I3.12p can be represented in the form 

dY Y^ R'^ / \ ~ 

— ^FyY- FrR + Fyy— + FyrYR + Frr— + O , ' , 

^ = LyY- LrR + Lyy^+ LyrYR + Lrr^ + O {\Y\^ , \R\^) ^ , (3.14) 

where the coefficients of the quadratic terms are the second-order derivatives 
with respect to the corresponding variables at the point of equilibrium Yq, Rq, 
and Fy=Ly-Sr> 0, Fr = Ly + Sr> 0. 

The matrix of the linear part of (|3.14p at the point of equilibrium has the 
following representation: 

\^ Ly —l^R J 
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with the characteristic polynomial 

\^ -trA- X + detA^O, (3.15) 

where ti A — Fy — Lji is the trace of the matrix A, and det A ~ FuLy ~ FyLji 
is the determinant of the matrix A. 

In the case when tr A < and det ^ > 0, we can argue that the system 
(I3.14P is stable in the linear approximation. We shall consider in more detail 
the situation in a small neighborhood of the boundary of the region of linear 
stability trA = ^, where ^ is a small sign-alternating quantity. This means that 
the quantities Fy and are sufficiently close to each other, and when they are 
exactly equal to each other, the divergence of the vector field of the variables 
Y (t), R{t) passes through zero. 

As det ^ > 0, we can set cj^ (/i) = dot A = FfjLy — FyL^ and Fy = Lj^ + fi. 

Equation p.l4p then takes the following form: 

- fiX + ul- fiLR^O, (3.16) 

Differentiating (|3.16|) with respect to the parameter /z, at /i = 0, we obtain: 

— = - - (3 17) 

d/i 2 2luq 

From (|3.16p . it follows that the real part of the derivative of the eigenvalue 
does not vanish. 

In other words, the eigenvalues cross the imaginary axis with a non-zero 
velocity. 

Thus, the conditions of Hopf 's bifurcation theorem are fulfilled, and, in the 
system (|3.14p , when the stability of the complex- focus-type state of equilibrium 
changes, the formation (or annihilation) of a bifurcation of the limit cycle takes 
place, accompanied by the generation of a corresponding self-oscillation regime. 

As a bifurcation parameter, we have chosen Fy — Lr, which is equivalent 
to the condition 

ly — By — Lji = 0. 

To analyze this bifurcation, we construct, using the change of variables Y = 
FrXi, R — LrXi + loqX2, the normal form of the system (|3.14p at /i = 0. 

As a result of transformations with the change of the time scale t = ujt, we 
obtain the system 

ii = -X2 + a20 Y + aiia;i2^2 + ^02-^' 

±2 = Xi + 620 Y + ^112^12^2 + &02Y, (3.18) 

where 

_ Fl^Fyy + 2FRLRFyR -f L^Frr 

020 — ; 
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On = pyr H ; ao2 — — ; 

i'R ^R 

^ _ Fr {FrLyy — Fyy) + '^FfjLjj (FrLyr — Fyr) + L\ {FrLrji — Frpi) ^ 

^{)Fr 

F'^Lyr — FrFyr + Ln {FuLrr — Frr) Frr 

Oil = ; 002 — ^RR — ■ 

^oi'R. i'R 

The system of two ordinary differential equations ()3.f 8p is the normal Poincare 
form; it can be directly applied for the evaluation of the main characteristics of 
the forming limit cycles, such as the amplitude, the frequency, and the period 
of oscillations; it can also be applied for the determination of the stability of 
periodic solutions. The relevant formulas are given work by V. Zang |19j : it is 
shown therein that the limit cycle can be either stable or unstable, depending 
on the values of the typical parameters. 

It would be in order here to emphasize that the above-mentioned work deals 
with the situation when the forming limit cycle is unique and has a quite definite 
type of stability. However, at the same time, the most important question in the 
studies of the Hopf bifurcation concerns the maximum number of limit cycles 
that can be generated from the state of equilibrium (the fixed point) under 
parametric excitation of the given system. This problem is completely resolved 
only for the quadratic case of polynomial systems by N. N. Bautin 6 : it is 
shown that the maximum number of limit cycles that can be generated in the 
quadratic system from a focus-type singular point is equal to three. 

To determine maximum multiplicity of the limit cycle in the system of dif- 
ferential equations (j3.18p . it is necessary to evaluate the first three Lyapunov 
focus quantities. 

Let us represent the system p.lSp in the form of a single complex differential 
equation in the variable Z — xi + 1x2, for 7^ 0: 

Z"^ - Z"^ 

Z = {i + fi)z + 920-^ +giiZZ + 320-^, (3.19) 

where 

9jk = gjk iajk,bjk) , j,k = 0,2, j + k = 2. 

The singular point turns from a focus into a center under the following 
conditions: 

1) 311 = 0; 

2) M = 320+311 = 0; 

3) = Im (320.911) = Im (3?i302) = Im (gfo3o2) = 0; 

4) ^ = 520 - 4gii = I502I - 2 |.9ii| = 0. (3.20) 

Expressions (|3.20p constitute the conditions for the existence of the first 
integral, or the Hamiltonian, of the system (|3.18p . In such a system, there 
exists an infinite set of periodic trajectories that continuously depend on the 
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initial conditions. Quite naturally, isolated closed trajectories (limit cycles) 
cannot exist under the Hamiltonian conditions. 

In the paper by H. Zoladek [48] , the following formulas for the evaluation of 
the three Lyapunov quantities for Eq. (|3.19|) are given: 

h = -^Im (520511) , 
^2 = - Im ((320 - 4gii) (520 + 511) 3iiff02) , (3.21) 

'3 = "^I"^ ~ l502|^) (320 +511 )5n 502 

Thus, using (|3.21[) . it is not difficult to establish the cyclicity of the singular 
point: 

1) no cycles are present if 

/i ^ 0; (3.22) 

2) a single cycle is present if 

// = 0, Im(52o5ii)^0; (3.23) 

3) two limit cycle coexist if 

/i = Im (520511) = 0, 520 7^ 4511; (3.24) 
3) three limit cycles are observed if 

M = 520 - 4511 = 0. (3.25) 
By virtue of the relation between p.lSp and p.l9p . i.e., 

520 = Y (020 - ao2 + 2611 + i (620 - ^02 - 2aii)) , 

511 = ^ (020 + ao2 - « (^20 + &02)) , 

we obtain a parametric restriction on the existence of three limit cycles. 
According to the condition (|3.25p . 

020 - ao2 + 2611 = 4 (a2o + ao2) , 

^20 - ^02 - 2aii = -4 (620 + &02) , 

or 

2611 = 3a2o + 5ao2, 

2aii =5620 + 3602. (3.26) 

Relations p.26p constitute algebraic conditions of the existence of three limit 
cycles in a system of the general form (|3.18p . They have been obtained by 
direct evaluation of the corresponding Lyapunov quantities for the six-parameter 
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Figure 3.2: Relative disposition of bifurcation surfaces in the neighborhood of 
the bifurcation I2 — 0. 

quadratic system without the use of any canonical models of the type of the 
system of Bautin, Andronova et al. [5], which employs a five-parameter form of 
the representation of considered models. 

In the presence of parametric perturbations, we write down the following 
differential equation characterizing the dynamics of the amplitude of oscillations 
with three small parameters Pj, j — 1, 3: 

p = /3lP + /32p2+^gp3^;^^4_ (3 27) 

Here, (3^ = fi, f]^ = h, = h, and ^3 ^ 0. 

The bifurcation diagram of the system for the case Z3 < is presented in 
Fig. 3.2 [4J. 

The plane N corresponds to bifurcation formation of a limit cycle from the 
fixed point O (the state of equilibrium). On the half-plane iV_ the loss of 
stability of the focus is " soft" , whereas it is " rigid" on N+ . The curved surface 
K corresponds to a non-local bifurcation of co-dimension one " double cycle" . 
The part Ki of this surface corresponds to a stable, from the outside, multiple 
cycle. The other part, K2, corresponds to an unstable, from the outside, multiple 
cycle. On the surface there is a rib of return C, i.e., a common line for the 
above-mentioned parts of the surface of the multiple cycle. On the line C, there 
forms a non-local bifurcation of co-dimension two, i.e., a bifurcation of a merger 
of three cycles. The line of the intersection of the half-plane and the surface 
Ki corresponds to a bifurcation of co-dimension "one plus one" accompanied 
simultaneously by a change of the stability of the focus and a merger of a remote 
pair of cycles. 

In Fig. 3.3, we present the bifurcation diagram in the neighborhood of 
(zero) for < 0. 

Returning to the description of the initial Keynes model (j3.ip . we want 
to emphasize that we have used a set of restrictions on the parameters of the 
linear part of the system (|3.14p that has allowed us to satisfy the requirements of 
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Figure 3.3: The bifurcation diagram for the bifurcation h = with I3 < 0. 

Hopf 's theorem on the existence of the hmit cycle. At the same time, we have not 
stipulated any restrictive conditions on the coefficients of the quadratic terms, 
and, therefore, we may argue that, in the most general case, the maximum 
possible number of limit cycles in the neighborhood of the state of equilibrium 
Yq, i?o is equal to three. 

Up to this point, we have studied the behavior of the system (|3.14p on 
one boundary of the stability region determined by the trace of the matrices 
of dynamics. In what follows, we shall assume that both the trace and the 
determinant of this matrix are small sign-alternating quantities. Let det A = 
—jjL-^ and tr A~ 

In this case, the characteristic equation (|3.15p takes the form 

- ^l2\ - Ml = 0. (3.28) 

For Ml = = 0, the eigenvalues are multiple and equal to zero: X1.2 = 0. 
Such degeneracy in the linear part of ()3.14|) may lead to the formation of the 
Bogdanov-Takens bifurcation [7]; a study of this bifurcation is the subject of 
what follows. 

As is obvious, in this case, one should choose two bifurcation parameter 
determined by the expressions for the trace and the determinant of the matrix 
A: 

Fy - Lfi = M2, 

FRLy-FyLR^-fi^. (3.29) 

Let the functions, found from condition (|3.29p . serve as the bifurcation pa- 
rameters: 

Fy = Lr + M2, 
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L'r + i_R/i2 - Ml 



(3.30) 



The next important stage is the construction of a corresponding normal 
form for the considered bifurcation of co-dimension two. Making the change of 
variables Y = Lnyi + 1/2, R = Lyyi, we obtain the following system: 



2/1 = ?/2 + m2oyl + "^112/1^2 + "^022/2- 



2/2 = fiiVl + M22/2 + n-202/1 + "112/12/2 + "022/2: 



where 



W20 



L\Lyy , r t , ^yLyR 
h LrLyr H 



mil 



2iy 

LrLyr 
Ly 



+ Lyr, mo2 



2 ' 
Ly Lyy 



"20 = I I -Fyy - -j—L 



JYY 



2LyLr \ Fyr — J— Lyr 



+Ly I Frr — -p-i 



RR. 



nil — Lr [ Fyy — j^Lrr^ 



Ly ( Fyr. — -j^Lyr. 



"02 — — \^RR — J~^^RR 

Carrying out a nonlinear reduction of the variables yi and y2, i.e., 

mil + "02 2 , 

yi=Ul-\ Ml + mo2UiU2, 



(3.31) 



2/2 = M2 - mo2Mi + no2MlM2, 

and dropping terms of order higher than two, we arrive at the following system: 



Ui = Ui, 

,2 



U2 — Ml^l + ^^2^2 + "20^1 + ("11 + 2m2o) UiM2- 



(3.32) 



By means of the substitution wi — ui + b, W2 — U2, we eliminate from (|3.32p 
the linear, in the variable U2, term: 



where 



Wi = W2, 

W2 ^ dl + O2W1 + 7l20W^ + (rill + 2m2o) WlU'2. 

"20 ..2 Ml/^2 



(3.33) 



(nil + 2m2o) 



2A*2 



nil + 2m2o ' 
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Figure 3.4: The diagram of the Bogdanov-Takcns bifurcation in the system 
(3.34) for 5 = -1. 



2n 



20 



-M2- 



nil + 2m2o 

To complete the construction of the normal form of the system p.3ip - p.33p . 
we need yet another scahng of the variables, ^i = ^ , ^i = "^"^^ W2, and 
of time, t^lKlT, K = "ii+2'"2o ^ q. 



where 



^1 ~ '>2i 
^2=^1+ a2Ci + Cl + '5''ClC2 



ai = (^^ll +2m2o) A'^6li, a2 = K^02, S" = sign (A') = ±1. 



(3.34) 



Eliminating 9i and O2, we express the small parameters Oi, a2 in terms of 
Ml, M2: 

a2 = ifVi - 2iv:/^2- (3-35) 

The bifurcation diagram for the case S* = — 1 is represented in Fig. 3.4. 

Our analysis of the system p.34|) we begin with the evaluations of the coor- 
dinates of the fixed point. As is obvious, since ^2 = 0, they are positioned on 
the horizontal axis in the phase plane and satisfy the quadratic equation 



ai +^2^1 +e? = 0. 



(3.36) 
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Equation (j3.36p may have from zero to two roots. The discriminant parabola, 

T={(ai,a2):4ai-a2 = 0}, (3.37) 

is related to the "fold" bifurcation. Along this line, the system p. 341) has 
equilibrium with the zero eigenvalue. If 02 7^ 0, the "fold" bifurcation is non- 
degenerate, and, on crossing the line T from left to right, two states of equilib- 
rium are formed. 

Denoting the left and the right states of equilibrium by Ei and E2 , respec- 
tively, we get: 




The point a = separates the two branches, T_ and T_|_, of the "fold"- 
bifurcation line for a2 < and 02 > 0, respectively. Along the line T_, a stable 
node El coexists with a saddle point E2, and, in the vicinity of r_|_, an unstable 
node El coexists with a saddle E2. 

The vertical line ai = is the line on which the state of equilibrium Ei has 
a pair of eigenvalues of the zero sum: Ai + A2 = 0. The lower part 

H = {(ai, 02) : ai = 0, a2 < 0} (3.38) 

is related to a non-degenerate Hopf bifurcation, whereas the upper half-line 
is a non-bifurcation line related to a neutral saddle. As a result of the Hopf 
bifurcation, a stable limit cycle is generated, i.e., the first Lyapunov quantity is 
negative. 

The cycle exists in the neighborhood of H, for ai < 0. 

The state of equilibrium E2 is still a saddle for all the values of the parameters 
to the left of the line T. Here, there are no other local bifurcations. 

Let us pass around the point ai = 02 = in circle of a small radius counter- 
clockwise. In region 1, both states of equilibrium and cycles are absent. When 
passing from region 1 to region 2 through the part T_ of the multiplicity line, 
two states of equilibrium are formed in the system (|3.34p : namely, a saddle E2 
and a stable node Ei. Further, the node turns into a focus that loses stability as 
a result of a Hopf bifurcation, which is accompanied by the formation of a stable 
limit cycle on the neutrality line of the focus H. This limit cycle disappears on 
the hue P = {(ai, q;2) : ai — —^a\, a2 < O}, being destroyed as a result of a 
global bifurcation on the loop of the separatrix of the saddle. (The period of 
motion along the cycle grows to infinity at that.) Finally, when passing from 
region 4 to region 1 through the part of the multiplicity line, the unstable 
node merges with the saddle, and both of them disappear. 

Using formulas (|3.35p . we write down the equations for the bifurcation lines 
in terms of the initial parameters fii, fi2- 

a) the line of the "fold" bifurcation is given by 



T = {(Mi,M2) : Ml = 0, H2 + 0}; 
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b) the line of the Hopf bifurcation is given by 

c) the hne of global bifurcation is given by 

r K 

P = { (Mi, M2) : (J^^2 - ^A^i) (7^2 - 6XAii) = 0, /i2 > yMi 

This point completes our study of the Bogdanov-Takens bifurcation of co- 
dimension two, as applied to the system p.l4p . The case S" = +1 can be 
studied using the substitution t —t, £,2 ~^ ~^2- Here, a substantial difference 
lies in the instability of the limit cycle. Besides, it is important to note that the 
limit cycle in the vicinity of the point of the Bogdanov-Takens bifurcation has a 
frequency which is proportional to the square root of the small parameter. This 
means that a business cycle in the Keynes model may have a very large period, 
which rather difficult to detect in numerical simulations of the initial model. 



3.3 Bifurcations in the nonlinear Kaldor model 

The main assumption of this model, constructed in the Keynesian spirit, is that 
the investment and savings functions are substantially nonlinear functions of the 
income Y and of the fixed assets rate K [JT] . As regards the investment function 
/ = I(Y,K), it is assumed that the limit propensity to invest ly = §y 
positive, although it is variable. This means that it takes the so-called "normal" 
value for "normal" values of the income rate Y. For values of income that are 
lower than a given "normal" interval, the limit propensity to invest declines 
as a result of losses of income in the period of low activity rate compared to 
the "normal" rate. It also decreases for values of Y that are higher than the 
"normal" interval because of a positive effect of the scale of expenditure and 
an increase in it. Thus, the investment function is an S-shaped curve. Besides, 
Kaldor assumes that a higher capital assets rate leads to a decrease in the 
marginal efficiency of the fixed assets; that is, ^ = —Ir < {Ik > 0). 

The savings function is also nonlinear: S = S {Y, K). The limit propensity to 
save By = ^ is positive and less than unity, although it varies. This assumption 
may be justified as follows: there exists a "normal" rate of the propensity to 
save that corresponds to a "normal" interval of changes in the income. Below 
this interval, the savings decrease towards consumption, whereas above this 
interval, they increase. In other words, the savings function is an upturned 
S-shaped curve. 

Additionally, Kaldor assumes that Sk = ^ is a positive quantity; that 
is, the savings function is shifted upwards with an increase in the capital rate. 
This assumption is questioned by a number of researchers [101 [Sj, because 
Kaldor himself did not provide any satisfactory justification for it. In our further 
consideration, we shall assume that Sk may change sign. 
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Figure 3.5: The structure of the states of equUibrium in the Kaldor model. 

In Fig. 3.5, we present quahtative dependencies of the investment and sav- 
ings functions on the income rate, for a fixed value of the capital rate. 
The dynamic Kaldor model is represented by the following equations: 



where Iq {Y, K) is the realized investment, which, generally speaking, is different 
from the planned one I ( Y, K) ; the quantity a > characterizes the rate of the 
change of the income in the time domain. 

The functions / (Y, K) and S (F, K) are symmetric with respect to the point 
of equilibrium. (Such an assumption is completely justified.) We also assume 
that the dependence of the above-mentioned functions on the variable / is 
linear, and that the realized investment coincides with the planned one, i.e., 
/q = /. It is then convenient to introduce new variables Y = Y — Yq and 
K — K — Kq that constitute deviations of the initial investment and capital 
rates from their equilibrium values at the point C. 

By symmetry, the investment and income functions are odd, and they can 
be represented in the following form: 



Y = a[{I {Y,K)~S{Y,K))], 



K = IoiY,K), 



(3.39) 





S (y, = SkK + SyY + SsY^ + O (f^) . (3.40) 
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Here, I3 and S3 are the corresponding coefficients of the Taylor series expan- 
sion of the initial function, and, besides, the condition ly > Sy holds. Thus, 
the system (|3.39p has the following explicit form: 

Y = a{Iy-SY)Y-a{lK-SK)K + a{h-S3)Y\ 

K ^ lyY ~ IkK + hY^ (3.41) 
The system p.4ip has three states of equilibrium: 

a) the trivial one, Yq = 0, Kq = 0, which corresponds to the point C in Fig. 

3.5; 

b) nontrivial ones, - T^^^g^, Ki,2 = %i^^§^Y,.,, which 
corresponds to the points A and B in Fig. 3.5. 

The points A, B and C represent possible variants of static equilibrium. 
In [IT], it is argued that equilibrium at the point C is unstable, whereas it is 
stable at the points A and B. At the point C, the instability of equilibrium 
is due to the fact that, for yA < y < Uc, savings exceed investments and a 
surplus appears in the goods market, which provokes a further decline in the 
production. In the case yc < y < ys, since the volume of investments exceeds 
that of savings, a deficiency of goods occurs, which stimulates a growth in the 
production. 

As regards the stability of the points A and B, we point out that a devi- 
ation from A or B to the right leads to a goods excess and to a decrease in 
their production, whereas a deviation to the left leads to a deficiency and to a 
growth in the production. The state of the economic system that corresponds 
to the point A is characterized by a low volume of investments that is insuffi- 
cient for a complete reimbursement of the worn-out capital. A decrease in the 
capital, after a certain time, will raise the entrepreneurs' propensity to invest, 
and the demand for investments will grow, which will lead to an increase in the 
investment function / {Y, K) ; equilibrium will be destroyed. 

On the contrary, the point B represents a state of equilibrium with high 
economic activities. As a result of an achieved optimal capital volume, the 
demand for investments starts to fall, the value of the function / {Y, K) begins 
to decrease, and the economy leaves the state of equilibrium. 

In the process of changes in market conditions, when the graphs of the 
savings and investment functions move towards each other, the points A and 
C may merge. In the opposite case, the points B and C may merge. It is 
important that these states of equilibrium, i.e., A, B and C, become unstable. 

In this regard, it is reasonable to study in detail the situation, when all the 
three states of equilibrium, A, B and C, are sufficiently close to each other and 
may merge into a single point. Our further consideration will be devoted to a 
detailed analysis of properties of the dynamic system (j3.4ip , taking into account 
the above-mentioned assumption. 

Consider the behavior of the system (|3.4ip in the vicinity of the trivial state 
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of equilibrium. The linear part of (|3.4ip has the characteristic polynomial 

+ {lK-a {Iy - Sy)) X + a [IySk + IkSy) - 0. (3.42) 

Let us assume that the coefficients in p. 421) are small quantities, i.e., = 
—a{lY Sk + IrSy) and ^2 = a (/y — Sy) — Ik- In this case, the system 
l|3.4ip may have a Bogdanov-Takens bifurcation, i.e., a so-called "double-zero" 
bifurcation I46j. As bifurcation parameters, we choose the following ones: 
iK^a [Iy - Sy) - and Sk = {W - Sy) - ^ + A^a- 

As fii and are small quantities, it is obvious that Sk < 0. Making the 
change of variables Y — [a {Iy — Sy) — ^^2] + ^2 and K — IyUi, we reduce 
the system to the form of a nonlinear oscillator. 

We have: 

Ui^U2 + aaoUl + aaiC/f J72 + ai2C/iC/| + aosUi, 
U2 = ^ilUl + Ai2C^2 + bsoU! + b2iUlU2 + h2UiUi + bosUi, (3.43) 



where 



0-30 = 7^ (a (Iy ~ Sy) - A*2)^ ; 021 = 7^ {a {Iy - Sy) - ^2)^ 

7y Iy 

ai2 = ^{a {Iy - Sy) - ^2) ; ^30 = 7^; 

lY J-Y 



bi] = 



a ( S'y - ^ ) + Ai2 



* + j = 3; j = 0, 3. 



To construct the normal Poincare form of the system of ordinary differential 
equations p.43p . we perform a nonlinear reduction of the variables of the state: 

U2^V2- azoVl + ^-fvlV2 + b^^ViVl (3.44) 

As a result of transformations, taking into account expressions p.44p . we 
obtain the following representation of the system (|3.43p : 

Vi = V2, 

V2 = fiiVi + Ai2^2 + MV{^ + NV^V2, (3.45) 

where M = 639, = 621 + Sasp. 

For ^1 = 0, /i2 = 0, we have: N = Sa^ {Iy - S'y)^ ~ S3 + Iy ~ Sy) , 

M = a"^ {Iy - Sy)^ (^^j^ ~ Ss'j . It is reasonable to reduce the system (P^iSl) . 
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Figure 3.6: The bifurcation diagram for 5* = +1. 

which is already a normal Poincare form, to a still simpler form by means of 
the linear transformation 

xi=p^/\M\Vi, X2=v'VW\y2. r=-t, 

P 

where 



N 


3 


1 1 


M 


a 


Iy — Sy Ssy — Syj^ 



Finally, we get: 

±1 = X2, 

±2 — p'^fJ-iXi + pfj.2X2 + Sxl — xlx2, (3.46) 

where S — sign {p) = ± 1 . 

As regards the system (13.461) . it is not difficult to notice that it is invariant 
with respect to the substitution xi — + ~xi and X2 — > —X2, and it always has 
trivial equilibrium Eq = (0,0). The two other possible states of equilibrium 
have the coordinates Ei_2 — (±-\/— O); they exist for fi^ < 0, if S ^ 1, and 
for > 0, if S' = -1. 

Important, for the system (|3.46p . is the fact that all the three states of 
equilibrium merge into the single trivial one for /ij^ = 0. 

Let S = 1. The bifurcation diagram is presented in Fig. 3.6. 

In region 1, there is the single trivial state of equilibrium Eq that is a sad- 
dle point. On crossing the lower branch of the line Ri — {(/^i,M2) ■ Mi — 0}, 
a "pitchfork" takes place, accompanied by the appearance of a pair of sym- 
metric saddles -Ei,2, until the trivial equilibrium Eq becomes a stable node. 
This node turns into a focus in region 2 that, on crossing the half- line i?2 ~ 
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Figure 3.7: The bifurcation diagram for S — —1. 

{{fJ-i, ^J■2) : Ml = 0, < 0}, undergoes a Hopf bifurcation accompanied by the 
generation of a stable hmit cycle. 

On crossing the line R3 = {(^1,^2) ■ M2 = "f^ij Mi < 0}> ^ global hetero- 
clinic bifurcation, accompanied by the appearance of corresponding orbits that 
are related to the saddles -Ei,2, takes place, and, in region 4, a heteroclinic cycle 
is formed. Further, all the three states of equilibrium coexist up to the crossing 
of the upper part of the straight line i?i, and a return to region 1 takes place. 

Consider the case = — 1. The corresponding bifurcation diagram is de- 
picted in Fig. 3.7. 

In region 1, there is the single trivial equilibrium Eq that is a stable node; fur- 
ther, it goes over to a focus. On the half-hne Li = {(M1JM2) • M2 = Mi < 0}, 
a Hopf bifurcation takes place, and a stable limit cycle is generated. Two un- 
stable nodes separate from the trivial equilibrium on crossing the upper part 
of the line L2 = {(Mi;M2) • Mi = 0> M2 > 0}' when passing from region 2 re- 
gion 3, as a result of a "pitchfork" bifurcation. In region 3, all the three 
states of equilibrium are localized inside a "large" limit cycle. On the half- 
line L3 — {(M1JM2) • M2 =PMi> Ml > 0}, the nontrivial focuses Ei^2 undergo 
a Hopf bifurcation. This bifurcation leads to the appearance of two "small" 
unstable limit cycles around the nontrivial states of equilibrium. The points 
of equilibrium become stable. Thus, in region 4, there are three limit cy- 
cles: an external "large" one and two internal "small" ones. Along the line 
L4 — {(M17M2) ■ M2 = ^/^i' (^1 ^ "small" cycles form a symmetric fig- 

ure that resembles the Bernoulli lemniscate, with the center at £^0i as a result 
of the occurrence of a global homoclinic bifurcation. Along the line ^4, the 
saddle Eq has two homoclinic orbits. These orbits can be transformed from 
one into the other by means of symmetry transformations. On crossing the 
line Li, when passing from region 4 to region 5, not only the "small" cycles 
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are destroyed, but also an external "large" unstable cycle is generated. Then, 
in region 5, two "large" cycles coexist: the external one is stable, while the 
internal one is unstable. These two cycles merge and disappear along the line 
L5 = {(Mi, /i2) ■ M2 — kop^J-i, Hi > 0}, where fcg — 0.752 . . .. This is a saddle- 
node bifurcation of the limit cycle. After the occurrence of this bifurcation, 
there are no limit cycles in the system. In region 6, all the three states of 
equilibrium are present: the trivial saddle Eq and two stable nontrivial focuses 
(nodes) £'1,2- The nontrivial states of equilibrium merge with the trivial one 
on the lower part of the line L2 = {(Mi, fJ-2) '■ f^i — 0, /I2 > 0} as a result of a 
"pitchfork" bifurcation, and we return to region 1. 

Thus, the behavioral properties of the system (|3.46p are determined by 
comparison between the quantities ly — Sy and ^3 — ^f^- If the condition 
ly — Sy > S3 ~ is fulfilled, then S — sign (p) = — 1. In the opposite case 
ly — Sy > S3 — ^j^, we get: S* = +1. Returning to the question of the choice 
of the sign of the quantity Sk , we note that states of equilibrium of the saddle 
type are possible in the initial dynamic system (j3.4ip only under the condition 
Sk < -^TT^, i.e., for Sk < 0. 

Thus, in this chapter, we have demonstrated dynamic behavior of the Kaldor 
model in its whole variety, in the case when three states of equilibrium degen- 
erate into a single one. We have also followed the hierarchy of instabilities 
accompanied by a cascade of corresponding bifurcations ^IjJ . 



Chapter 4 



Dynamics of economic 
processes with a lag 

In many economic models, when constructing balance relations, one has to take 
into account a lag responsible for various interactions. In Chapter 2, we have 
already considered the Goodwin model involving a lag of two types: on the part 
of the demand for investments there is a lag with fixed duration of the action 
of the accelerator, whereas on the part of the supply there is a continuously 
distributed lag. In this chapter, we shall demonstrate the application of methods 
of the theory of nonlinear oscillations in the class of differential equations with a 
lag and of integro-differential equations that may have a continuously distributed 
lag. 

4.1 Instability of price dynamics in Fisher's model 

In modern economic literature, there is a sufficiently detailed qualitative de- 
scription of the mechanism of the formation of domestic prices of the output 
based on an analysis of the dynamics of financial flows in export-import opera- 
tions. At the same time, special attention is paid to such an important factor as 
an index of the trade balance whose surplus or deficit determines the direction 
of price change. 

Before proceeding with the consideration of a formalized mathematical model 
based on classical Fisher's macroeconomic equation, it is necessary to put for- 
ward a number of assumptions [41j : 

1) a free-trade scheme, without any infiuence of governments and monopo- 
listic structures, is considered; 

2) the rate of the national income is considered to be sufficient, and the price 
rate is determined on the basis of qualitative theory of the currency of money; 

3) in the course of the considered period, changes in the supply of money 
stock are stipulated exclusively by a surplus (or deficit) of the balance of trade; 
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4) since exchange rates are considered to be fixed, they can be set equal to 
unity, which is equivalent to complete unification of international operations; 

5) transportation expenses, insurances and other costs are not taken into 
account either for commodity or financial fiows. 

In what follows, we use the following notation: Q is money supply; V is the 
rate of the currency of money; Y is the national income rate (a constant); P is 
the index of the domestic price; Pm is the index of the export price (a constant); 
M is the volume of import; X is the volume of export. 

The basic equation of Fisher's model has the form 

QV ^ PY. (4.1) 

The volume of export is a decreasing function of the domestic price: 

r)X 

X^X{P), — <0. 

On the contrary, the volume of import is an increasing function of the domestic 
price: 

dM 

M = M{P), — >0. 

A relation of the form P*X (P*) — PmM (P*) = characterizes a condition 
of trade balance equilibrium, and it is assumed that this algebraic equation has 
positive solutions that determine equilibrium values of the domestic price P* . 

A disturbance of equilibrium is accompanied by changes in the money supply 
and is expressed by means of the following equation: 

^^PX{P)-PmM{P). (4.2) 
at 

From expression (14. ip . it follows that a change in the money supply leads 
to a change in the domestic price P = P{t). We assume that this change is 
not instantaneous, i.e., there is a time lag determined by a constant positive 
quantity r. In this case, equation (|4.ip takes the form 

dP{t) _V dQ{t-T) 
dt Y dt ' ^ ' ' 

From (|4.2p and (|4.3p . we obtain a differential-difference equation that de- 
scribes the dynamics of the domestic price: 

^ = ^{P{t-T)X[P{t- T)] - PmM [P {t - r)]} . (4.4) 

Let us assume that the export function X (P) and import function X (P) 
are inherently nonlinear and can be expanded in a Taylor series up to the third 
power in the neighborhood of the state of equilibrium P*, i.e., 

X (P) = + (P - P*) + ^2 (P - P*f + (P - P*f + O (P^) , 
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M (P) = Mo + Ml (P - P*) + M2 (P - P*Y + M3 (P - P*y + O (p4) , (4.5) 

where Xi = ^ ^Qpi ^ , Mi = ^Tfgpr^: with i = 0, 3, are corresponding derivatives 
of the functions X (P) and M (P) at the point of equihbrium P* . We introduce 
a new variable P (t) = P (t) — P* that has the meaning of a deviation of the 
domestic price from its equihbrium value. In this case, taking into account (j4.5p . 
equation (|4.4p reduces to the following: 

^ = ^ [GiP(f- 1) + G2P2 ^ g^p3 1) ^ (5 (p4)] ^ (4 g) 

where i = rt, d = + P*X, - PnMi, i = T73. 

It is reasonable to begin an analysis of the process described by (|4.6p with a 
study of the conditions of local stability, restricting ourselves only to the linear 
part, i.e., 

,4.7, 

The characteristic equation for (|4.7p is given by 

A-^e-=0. (4.8) 

Using a well-known result of the theory of stability of differential-difference 
equations \43\ , as applied to ()4.7|) , we obtain the necessary and sufficient condi- 
tions of linear stability: 

0< ^<2- (4.9) 

From the form of the left-hand part of the double inequality (14. 9|) it follows 
that the quantity Gi is negative, whereas the right-hand part of (|4.9p sets the 
upper bound for the absolute value of Gi. 

Condition (j4.9p has a rather transparent and meaningful interpretation. To 
demonstrate this interpretation, we perform a transformation of the initial pa- 
rameters of the considered model (j4.6p . 

Let 

Gi = Xoil-rjx - Vm] , 

where 

_ P*Xi _ P*Mi 

Vx — ) Vm ~ —r? j 

Aq Mq 

under the condition P*Xo — PmMq. 

The quantities r]x and r|J^J are elasticities of the export and import functions 
with respect to the price P. As Xq > 0, the condition Gi < is equivalent to 
Vx + Vm > 1; which corresponds to the so-called Marshall-Lerner conditions 
[51]. At the same time, condition (|4.9p reduces to 

0<..+.M<l + ^. (4.10) 
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Thus, the economic interpretation of the conditions of local stability consists 
in the fact that the sum of the elasticities not only must be larger than unity, 
but it also must be smaller than an additional critical value. In other words, 
instability in the considered economic model may arise not only when the sum 
of the elasticities is sufficiently small, but also in the case when its value is much 
larger. 

Let us study behavioral properties of the initial dynamic system ()4.6|) in a 
small neighborhood of the bounds of the inequality (|4.10p . First, we consider 
the situation when stability is lost at the lower bound. To this end, we introduce 
a small parameter i>i = 1 — rjx — Vm- 

In this case, when the sign of vi changes, the eigenvalue of the linearized 
problem passes through zero, and a stationary value P* may either not exist or 
split into several stationary states. This means that a bifurcation of stationary 
solutions takes place. 

The differential-difference equation (|4.6p can be represented as follows: 



HP 

= A^P [t- 1) + A2P^ {i- 1) + ^3^' (t~ 1) , (4.11) 

where Ai — ^Gi, i — 1,3. 

One should bear in mind that the quantity Ai is small, i.e., Gi = XqVi. 
Additionally, we assume that A2 is also small if we introduce a small quantity 
U2 = Using the techniques of the central manifold method [35], one can 
prove that the differential-difference equation (|4.1ip . under the condition that 
Ai, A2 are small and the time lag is finite, is topologically equivalent to the 
differential equation 

HP 

= AiP (t) + A2P^t) + A^P^t) (4.12) 



in the neighborhood of P = 

By means of the lii 
represented as follows: 



By means of the linear change of variables P = P + , equation (|4.12p is 



dP 



= ai+a2P + a3P^, (4.13) 



where 



2^_^A^ _ 

27 AI 3A3 ' ' 3^3' 



The transformation P {t) — \j3\W {t) yields an explicit form of the normal 
Poincare form for the differential equation (|4.13p : 

'^^P^+P^W + SW\ (4.14) 



where 



/?-^^, Pi^rkv 1^2^012, ^ = sign/3 = ±1. 



/A. 
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Figure 4.1: The diagram of the " cusp-of-the-beak" bifurcation. 



For definiteness, we set 5 = —1. 

Equation (|4.14p can have three states of equihbrium. A " fold" bifurcation is 
determined by a curve R on the plane /Jj, given by a projection of the line 



A 



02 - ^ 



to the parameter plane. By eliminating V from these equations, we obtain the 
projection: 

i?={(/3i,/32) : 4/33 + 27/3? = 0}. 

The curve R is called a semicubical parabola, and it has two branches Ri, 
R2 that meet tangentially at the "cusp-of-the-beak" point (a cusp bifurcation), 
for Pi = 1^2 = 0. The corresponding bifurcation diagram is represented in Fig. 
4.1. 

In region 1, in front of the boundary line, there are three states of equilib- 
rium: two stable states and one unstable state. In region 2, behind the separa- 
tion line, there is a single state of equilibrium which is stable. A non-degenerate 
"fold" bifurcation takes place on crossing cither Ri or R2 at any point of the 
parameter plane Pi,(32y except for the origin. If the curve is crossed on passing 
from region 1 to region 2, the right stable state of equilibrium merges with the 
unstable one, and both of them disappear. Analogously, the left stable state of 
equilibrium merges with the unstable one on the line R2. 

On approaching the "cusp-of-the-beak" point, in front of the regionl, all 
the three states of equilibrium merge into a single one as a triple root of the 
right-hand side of the initial equation (|4.14[) . Of importance is also the fact 
that, in the course of the transition from a stable regime to an unstable one, the 
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phenomenon of hysteresis is observed in (|4.14p . and a catastrophe occurs [15] . 
The case S — I can be considered analogously. 

Let us now clarify the situation when the sum of the elasticities rjx + Vm is 
close to the right bound of the inequality (|4.10|) . We introduce into consideration 
a small parameter ^ ~ — (Ai + Then, the differential-difference equation 
(|4.6p can be transformed into 

dP 



^_ ^ii+-^P{t-l) + A2P'{t-l) + A:iP'{t-l). (4.15) 

The characteristic equation for the linear part of (|4.15p is given by 

A+(/i+|)e-^ = 0. (4.16) 

We should find out when this equation has a pair of purely imaginary roots 
A = ±iuj, = -1, w > 0. 

If A = iiw, then the conditions 

^/i + — ^ cos a; = 0, — ^/i + — ^ sintj — 

hold. 

From this fact, it follows that, for /i = 0, equation (|4.16p has a pair of purely 
imaginary roots for a; = ^. It is not difficult to show that (|4.16p has no roots 
with positive real parts. 

As A is analytic with respect to the differentiation of (|4.16p . at /i = 0, 
yields: 

d\ _ f +z 
9m ~ ^ + 1' 

Thus, all the conditions of Hopf's theorem on the existence of periodic solu- 
tions are satisfied, because the real part of the derivative of the eigenvalue with 
respect to the parameter does not vanish. 

Based on the above-mentioned considerations, we shall demonstrate that Eq. 
(|4.15p has a family of periodic solutions P^ (t) {e > 0) , where e is a measure of 
the amplitude maxf \P^ (f)|, and e is sufficiently small at that. 

The problem reduces to a study of the bifurcation of the generation (anni- 
hilation) of a cycle in the differential-difference equation (|4.15p . To reduce this 
functional equation to a complex differential equation, we use the method of the 
central manifold [37 . 

Equation (j4.15p contains a number of parameters; therefore, to simplify fur- 
ther consideration, we make the substitution P^ (t) = (t). 

For fj, = 0, equation (|4.15p takes the form 

^^-l{u{t-l) + u' (f- 1) +7^^ it- 1)) , (4.17) 

where 7 = 
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By the central manifold theorem, equation (|4.17p reduces to a differential 
equation in a complex variable: 



where 



• ITT - Z^Z 



521 = 27r 



520 = -511 = 502 = T^D, 

2 - Hi 3jTr\ - 7 - 

] D + -DD + iD^ 

5 4/3 



D = 



1 



D = 



1 



(4.18) 



(4.19) 



The existence of concrete values of the coefficients of the nonlinear part of 
Eq. (|4.18p allows us to use the formulas of [37] in order to determine stability, 
the direction of generation, the period, and the asymptotic form of periodic 
solutions of small amplitude of the limit cycle that realizes the Andronov-Hopf 
bifurcation from the stationary state. Using (|4.19p . we obtain an explicit form 
of the first Lyapunov quantity: 



Ci(0) 



1 



f 15 



2 TT /II 377r\ 

--2 T + — 



vr 

^'5 



11 



377r 



- + ^ . (4.20) 



The real part of ()4.20p is negative for 

16 - 447r 



7 > 7o 



157r2 



-0.826.. 



(4.21) 



This means that the limit cycle is stable if 7 > 79, and it is unstable if 
the condition (|4.2ip is not fulfilled. For the stable limit cycle, the following 
expressions for the major characteristics are derived: 

1) the amplitude is given by 



e = 



20/i 



1577r2 + 447r - 16 



2) the period is given by 



= 4 1 



3) the asymptotic form of the periodic solution is given by 



Ue {t) = 2e cos 



nt 



2e^ { — sin (nt) cos (vrt) — 1 

5 5 



At the same time, the cycle is generated in the direction /i > 0, and the 
emerging periodic solution is asymptotically stable. The corresponding regime 
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of the generation of self-oscillations is called soft. On the contrary, if the condi- 
tion 7 < 7g is realized, an unstable limit cycle takes place. The loss of stability 
with the generation of self-oscillations occurs rigidly, i.e., a sharp transition 
(jump) to a new stationary regime. In a realistic system, such a loss of stability 
results in a catastrophe. 

The most complicated behavior of the initial system (|4.17p is exhibited in 
the situation when the parameter 7 is close to its critical value 7q, i.e., when 
the quantity C = 7 ~ 7o is small. In this case, one can observe in the considered 
system the so-called Bautin bifurcation [37] that is characterized by a possibility 
of the coexistence of both the stable and unstable limit cycles. 

To analyze qualitative properties of the above-mentioned bifurcation, we 
should expand the right-hand side of (|4.17|1 in a Taylor series up to fifth-order 
terms. After that, we should employ the central manifold method to reduce 
the functional equation to a complex differential equation involving nonlinear 
fifth-order terms, which is necessary for the evaluation of the second Lyapunov 
quantity C2 (0) ^ 0. 

In this case, the complex differential equation has the form 

Z = Z {iiu + si + £2ZZ + CaZ^Z^) , 

where e\ — e\ £2 — £2 (0- 

Depending on the sign of £1, £2 and C2, the following scenarios are possible: 

1) C2 < 0, £2 < 0. When e\ passes from negative values to positive ones, 
the system softly achieves a stable self-oscillation regime; 

2) C2 < 0, £2 > 0. When £1 passes from negative values to positive ones, the 
system rigidly achieves a stable periodic self-oscillation regime. It is generated 
before the loss of stability by the state of equilibrium, together with an unstable 
oscillation regime that settles on the state of equilibrium at the very moment 
when stability is lost; 

3) C2 > 0, £2 < 0. The loss of stability is soft. However, the generated cycle 
is quickly annihilated in the process of merging with an unstable one, coming 
from a distance. After that, a new regime is rigidly excited in the system; 

4) C2 > 0, £2 > 0. This is classical rigid excitation. 

Consequently, whatever the sign of C2 (0), for corresponding sign of £2, our 
analysis reveals a qualitatively different, compared to the one-parameter case, 
phenomenon: for C2 (0) < there exists a rigidly excited stationary regime, 
whereas for C2 (0) > a softly excited regime turns out to be short-lived. In 
order to establish which of the two cases (C2 (0) < or C2 (0) > 0) is realized 
in reality, one has to perform scrupulous evaluation of the second Lyapunov 
quantity, which is in itself a rather good exercise in symbolic transformations. 

Thus, as a result of the study of dynamic properties of the differential- 
difference equation (|4.4p in a small neighborhood of the boundaries of the region 
of local stability, we can arrive at the conclusion that there exists a bifurcation 
of co-dimension two (this fact is in itself far from being trivial): 

1) on the left boundary (r\x + ?/m = 1), a " cusp-of-the-beak" bifurcation 
takes place; 
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2) on the right boundary (77 + rjj^j = 1 + Vr )' ^ Bautin bifurcation of 
the Umit cycle takes place. 

In conclusion, we want to point out that the application of mathematical 
methods to an analysis of concrete objects is associated with numerical results 
and a corresponding meaningful interpretation. In this sense, the role of qual- 
itative theory of differential-difference equations is somewhat different: it puts 
stress on a search for characteristic features of the phenomenon as a whole, on 
qualitative forecasts of its behavior. The objectives of the authors are the de- 
termination of irreducible topological structures that form the phase portrait of 
the system. The applied part consists in establishing a correspondence between 
these structures of the phase space and the considered economic processes, to- 
gether with carrying out a bifurcation analysis. At the same time, we have to 
take into account the properties of the realistic object that impose restrictions 
on both the phase variables and the parameters of the initial equations [9j . 

4.2 The cyclicity of innovation processes 

Presently, the humanity is concerned with a search for innovation ways of stable 
development of civilization. A new paradigm of the 21st century, i.e., the con- 
cept of stable development, has explicitly systematic, synergetic character. By 
stable development we should understand a synthesis of the necessities of stable 
economic, ecological and social evolution, which is realized simultaneously on 
global, national and regional levels; and, as is known, simultaneous cooperative 
action is the essence of the synergetic effect. 

In this work, by economic development we shall understand a substantially 
nonlinear process, characterized by spasmodic transitions from one stationary 
state to another. A fundamental basis of economic (innovation) forecasting is 
formed by prediction theory of N. D. Kondratcv [24 and innovation theory of 
Y. A. Shumpeter [47j . further developed by the modern Russian economist Yu. 
Yakovets [38] . 

The main theoretical prerequisites for the justification of qualitative fore- 
casting are the following: 

1) a prediction of economic innovations is based on accounting for the inter- 
action of the laws of statics (that determines a multiple balance of the function- 
ing of the economic system), of cyclic dynamics (a combination of the observed 
cycles of various duration), and of sociogenetics (the laws of heredity, variability 
and selection in the dynamics of technological and social-economic systems); 

2) depending on the period and the intensity of influence on the economy, 
evolutionary cycles differ substantially: there are medium-term, long-term and 
super-long-term cycles. Innovation oscillations facilitate technological crises ac- 
companied by changes in the structure of innovations. At the beginning of a 
super-long-term cycle, the most important seminal innovations are initiated; 
once during a few centuries, they radically change the structure of the economy 
by forming a new technological way of production. Kondratev's half-century 
cycles are characterized by basic innovations that determine the competitive- 
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ness of the production in the framework of a given technological formation of 
the economy. Medium-term cycles (of 10-12 years of duration) are caused by 
changes in the prevailing generations of technology; they are realized in the 
cluster of basic innovations and in the wave of improving ones. Thus, we may 
argue that innovation-technological cycles of various periods serve as the basis of 
economic cycles of corresponding duration: medium-term cycles (of 10-12 years 
of duration), long-term (Kondratev's) cycles, and super-long-term (civilization) 
cycles [39j . 

A rather detailed description of the phases of the innovation cycle and their 
application to an analysis of dynamic structure of the economy is contained in 
the monograph [26] . Conventionally, as the beginning of the cycle, we shall take 
a decrease in the efficiency of the prevailing generations of technics (technol- 
ogy) that leads to a decrease in the rate of economic growth and a drop in the 
standard of living of a substantial part of the population. This stimulates sci- 
entific and engineering activities in the direction of obtaining new technological 
solutions. 

However, conditions for their innovative application are generated only at 
the end of the depression phase and in the phase of revival, when a renewal of 
fixed capital takes place and the volume of investments grows, which stimulates 
a demand for innovations. The rate of economic growth accelerates, employment 
increases, and a growth in the standard of living of the population is observed. 
By the end of the revival phase, these processes reach their maximum. However, 
in the stability phase, against a background of a considerable income rate, the 
rate of growth falls until a crisis causes a sharp drop in the rate of growth 
and a decline in the gross domestic product (GDP), which substantially reduces 
investments. 

In the monograph by S. Yu. Glaziev [17j . it is pointed out that, simulta- 
neously with the acceleration of economic development, the influence of coun- 
teracting factors increases: this conclusion is drawn on the basis of multiple 
statistical measurement of the dynamics of GDP for different countries. As a 
result, economic growth either stabilizes or acquires cyclic character. These 
arguments are illustrated by a mathematical model that reflects feedback be- 
tween the rate of economic growth and the rate of the gross national product 
per capita. It is assumed that the influence of counteracting factors increases 
with cumulative growth of GDP. In other words, we assume that the rate of the 
growth of the cumulative volume of the production of a new product depends on 
the average weighted cumulative production volume in the past and not only on 
its volume at the moment. As a result, the dynamics of the production volume 
can be described by the integro-differential equation 



where yo (t) is the cumulative volume of the production of the new product; 

Q (t) is a function characterizing cumulative growth of the production; 

r > is a parameter that has the meaning of a technological limit of pro- 
duction growth. 




(4.22) 
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As a rule, as Q {t) , one employs a function with a corresponding set of 
characteristic time lags. In terms of the theory of automatic control, Q {t) 
is an impulsive transition function of a linear control system with an input 
signal Uo (t). The output of such a system is given by the quantity uq (t) — 
J-oo yo {t + s) Q (— s) ds, which is a convolution of the functions yo (t) and Q (t). 
In what follows, we shall use an operator representation of Q (t) in the form of 
a fractionally linear function of arbitrary order, namely, 

A +a„_iA +... + aiA + ao 
under the normalization condition 



f Q{-s)ds = l. 

J —OO 



The nonlinear integro-differential equation (j4.22p has two states of equiHbrium: 

1) vo = 0; 

2) Vo = r. 

The state of equilibrium j/q = is of no practical interest for our analysis, 
because it implies complete absence of the GDP. The second singular point, Uq = 
r, characterizing the limit of GDP growth, is a major economic factor; hence, 
our objective will be a study of dynamic properties of the process described 
by the integro-differential equation (|4.22p in the neighborhood of this state of 
equilibrium. 

Let us rewrite Eq. (j4.22p in terms of a new variable y — yo ^ r that has the 
meaning of a deviation of the production volume from its equilibrium value: 

/o 
y{t + s) Qn i-s) d~y y{t + s)Qn (-s) ds. (4.24) 

The characteristic equation of the linear part of (|4.22p has the form 

A + rg„(A)=0, (4.25) 

where Qn (A) is defined by expression (|4.23p . 

Let us consider the simplest case n — 1, Qi = j^^. To satisfy the nor- 
malization condition, it is necessary that feo,ao > 0. Then, equation (|4.25[) is 
represented in the form of the quadratic equation 

A^ + oqA + rba = 0. 

Note that, for Qq > 4r6o, the state of equilibrium y* = is a stable node. 
However, for Qq < 4r6o, in the vicinity of the equilibrium point, stable oscil- 
lations with the attenuation rate ao are observed. This type of equilibrium is 
called a stable focus. In this system, a transition from a node to a focus has no 
bifurcation character, which means that the system is stable. 
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In the sense of the diversity of dynamic behavioral properties, the situation 
with n = 2 is more interesting. In this case, 

^'^ ^ > ~ T2~^ T"^ — ' 

A +aiA + ao 

and the spectral equation (|4.25p takes the form of the cubic equation 

+ aiX^ + {ao + rbi)X + rbo = 0. (4.26) 
If the coefhcients of Eq. ()4.26|) satisfy the relation 

ai{ao + rbi)^rbo, (4.27) 

we obtain the solution Ai 2 = ±iuj, A3 — — ai, where lu^ — , , r — , , 

— 1. As the parameters uj and r are positive, for the coefficients we have: 

oi, oq, bo > 0, and bo > aibi at that. 

The presence of a pair of purely imaginary eigenvalues in the spectrum of 

the linear part of the integro-differential equation (|4.24p implies a possibility of 

the excitation of a self-oscillation regime, i.e., an occurrence of a limit cycle as 

a result of an Andronov-Hopf bifurcation. 

In order to transform the integro-differential equation (|4.24p into a system 

of nonlinear differential equations, we make the change of variables Xi = y, 

X2 = uo — r, 0:3 = Up. As a result, we arrive at the following system of ordinary 

third-order differential equations: 

±1 = -rx2 - X1X2, 

±2 = X3, (4.28) 
±3 = boxi - (flQ -I- bir) X2 - axxj, - biXiX2- 

Now, we shall demonstrate the application of Hopf 's bifurcation theorem to 
the system of autonomous differential equations (I4.28p . 

First of all, we should decide on the choice of the bifurcation parameter 
whose critical value allows for the existence of purely imaginary eigenvalues, 
i.e., of those whose real part vanishes (ReAi_2 = 0), whereas the imaginary 
part is nonzero (ImAi.2 = ±tj). It is convenient to take r as the bifurcation 
parameter and to study the properties of (|4.28p in a small neighborhood of the 
critical value = ba'^-aibi • That is to say, we consider the value r = + /i, 
where ^ is a small quantity. In this case, the eigenfrequency of the system (|4.28p 
in the linear approximation is defined by the expression 

^.^aobo^boro ^ ^^^9) 
Oq — ai&i ai 

In the new notation, the characteristic equation (|4.26p takes the form 



A^ + aiA^ -1- (w^ + bifi) A + feoM + ^lUj'^ = 0. 



(4.30) 
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For /i = 0, equation (|4.30p has the above-mentioned solution 

M,2 = ^i^; A3 = — ai (Re Ai^2 > Re A3) . 

Let us differentiate Eq. (|4.30p with respect to the parameter /i. For ^ — 0, 
we obtain: 

dX {bo — aibi) uj + i (^aibo + biLu'^) 



A'(0) 



Re A (0) = —-^5 ^, IniA (0) = - — ^. (4-31) 

Thus, all the conditions of Hopf's theorem are fulfilled, and we may argue 
that there exists a bifurcation of the generation of a cycle from a complex focus. 
Using the change of variables 

ai ^2 2/2 , 

xi = 7-2/1 + -T-ys, X2 = h ys, xs^yi- aiya, 

bo bo u 

we represent the system of nonlinear differential equations (|4.28p in the form 
that is convenient for the construction of the normal Poincare form. As a result 
of the transformation, we obtain (/i = 0): 

2/1 = -W2/2 +Fi (2/1,^2,2/3), 

2)2=^^2/1+^2(2/1,2/2,2/3), (4.32) 
ya = -aiys + -F3 (2/1, y2, ya) , 

where 

F^ (yi,y2,y3) = ^j9'(2/i,2/2,ya) , « = 1,3; 
, ^ ai ai u; 2 

<^(yi,y2,y3) = 7 — 2/12/2 + ^-yiys + 7-2/2y3 + -r-v?,; 

bouj bo bo bo 

^1 = -2wImA'(0), A2 = -2tjReA'(0), yl2 = -2 Re A' (0) . 

Let us reduce the order of the system of differential equations by means of 
the introduction of the new coordinates 

z = yi+iy2, z = yi-iy2, i^ = y3- 
As a result, equations (|4.32p are represented as follows: 

i = iujz + G {z, z, v) , 

v = —aiv + H {z,z,v) , (4.33) 
G (z, z, v) = Fi (z, z, v) + iF2 (z, z, ;/) , (z, z, ;/) = F3 (z, z, ;/) . 
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For the sake of a further analysis of the main characteristics of the hmit 
cycle, we use the central manifold method, which yields the relation 

2 -2 

v = W{z,z):^ W2.Q— + wxxzz + wq2^ + (|z|^) , (4.34) 

where 

W20 = (ai + 2iw)~"^ /i20, wii=a^^hii, wo2 = (ai - 2iw)""^ /io2, 

As a result of the substitution of expressions (I4.34p into (|4.33p . after certain 
necessary transformations, we obtain 



where 



z = iLoz + G20— + Giizz + Go2y + G12— + . . . , (4.35) 

2a; 

_ -fli (60 - ai&i) (t^ (2&0 - fli^i) + i {aibo + 2&iw^)) 
" 8cc;62 (a;2+a2) (4^2 + ^2) • 

Now we possess all the necessary data for the evaluation of the first Lyapunov 
quantity: 

Ci(0)-^ + ^|Go2|'. (4.36) 
From and it follows that 

Re C, (0) = Tj'"r1wfrt^^ < 0, (4.37) 
86§(w2+a?)(4u;2 + a?) 

since bo > aibi and, accordingly, 26o > ai^i. 

Let us consider peculiarities of the limit cycle that is generated from a com- 
plex focus when a pair of roots cross the imaginary axis. A stable focus takes 
place when fi < 0. When /i passes through zero, the focus at the origin loses its 
stability. For n — 0, the focus at the origin is stable but non-coarse: the phase 
curves approach zero exponentially. 

For /i > 0, the phase curves, having moved away from the focus at a distance 
proportional to fi^ , get wound round the stable limit cycle. In other words, 
the loss of stability under the change of the sign of /i is accompanied by the 
excitation of a stable limit cycle whose radius grows as /i^ . 

Thus, the stationary state loses stability, and a stable periodic regime is 
generated in the direction /i > 0; its amplitude is proportional to the square 
root of the deviation of the parameter from its critical value. Corresponding 
excitation of self-oscillations is called soft. 
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Concerning the obtained limit cycle, it is not difficult to derive explicit ex- 
pressions for its main characteristics following the methods of [37] : 

a) the amplitude is given by 

b) the period is given by 

27r 



T^-{l + r,e' + 0{s')); 



The periodic solution itself, up to the choice of the initial phase, takes the 
form 

yi = Re z, 1/2 = Im z, = Re (w2oz^) , 



Z = ee~ + 



(e-^ -3e^) +0(e^) . (4.40) 
From (|4.40|) . it follows: 



27rA ^iflie^ ( A_-Kt\ AaOie^ ( A_-Kt\ 



y2 it) = ssm ( — - ^^cos — - ^^sm — ) , (4.41) 



2/3 it) = -4t^ i-yl W + -yi (t) (0 + y' (*)) ■ 

Returning to the initial variables, we obtain an approximate form of the 
solution for the cumulative volume of the production of the new product: 

2 

2/0 it) = r + ^yi it) + i^yg it) ; (4.42) 

OO Of) 

for its relative speed (rate) of growth, we have: 

yoit) f,. y2 it) , . . . 

- -^2 (0 = ^ - 2/3 it) . (4.43) 

The above-described mechanism of the generation of an economic cycle has a 
number of characteristic features that cannot be explained within the framework 
of classical linear theory. 

Firstly, in contrast to the linear model, economic oscillation processes in the 
form (|4.4ip - (|4.43p do not possess symmetry. This means that, inside the cycle, 
its phases may be different both with respect to their form and meaning. In 
particular, expansion and decline have different duration. Secondly, according 
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to (|4.38p . the amplitude of oscillations depend on internal parameters of the 
model, and it is not an entirely external characteristic. 

Thirdly, the period of the cycle is a function of its amplitude, which, to a 
certain extent, explains the irregularity of cyclic dynamics. 

Furthermore, the evolution curve of the cumulative national product is a 
superposition of oscillations with different frequencies with respect to a certain 
trend. Such a trend can be represented, for example, by a logistic curve: in our 
case, logistics easily manifests itself when the lag Q {t) is switched off. 

As an illustration of the justification of the application of the mathematical 
models of economic development employed in this study, we should remind 
the reader of the statement of Y. A. Shumpeter quoted in [17]: "Since a long 
time ago, experts in political economy have a habit to mention 'trends' and 
'lags'; they are probably aware of the fact that businessmen react not only to 
given quantities but also to the rate of their changes, not only to the existing 
quantities but also to those that are expected in the future. In any case, in the 
recent years, exact theories of lagging adaptation, of expected actions, etc. have 
been developed. Technical means have been developed or borrowed from other 
fields. As regards the latter, the introduction to economic theory of functional, 
developed by Vito Volterra, has been the most important event... As it seems, 
the new methods point to the possibility of a colossal variety of wave-like motions 
in the economic life that can be applied to the explanation of cycles without 
any reference to the principle of the realization of new combinations... These 
new techniques of the analysis substantially extend our abilities to explain the 
forms of the manifestation of reality...". 

The results obtained in this study agree with the laws of cyclic dynamics: 
they demonstrate the irregularity of the evolution of the economic system, a 
periodic change of the phases of the cycles of complex systems and a change of 
the cycles themselves |21) . 

This, in turn, implies the relevance of the prediction of cycles and of cri- 
sis phenomena, of timely detection and identification of negative tendencies of 
development for the methodology of social-economic forecasting. At the same 
time, one should always bear in mind that cyclicity is a general property of the 
behavior of a large number of dynamic systems of various nature, and, when 
making extrapolation to the future of the tendencies formed in the past, it is 
necessary to be able to apply powerful modern techniques of the theory of non- 
linear oscillations. This approach allows us to provide a high quality analysis 
of crisis states of the studied objects and systems, it also facilitates a search 
for optimal, efhcient ways of leaving these states in order to create anti-crisis 
programs of economic and social development of the society. 



Instead of the conclusion 



Nowadays, the theory of economic cycles seems to be the most disputable section 
of macroeconomics. What kinds of economic or non-economic factors generate 
oscillation motion? What is the mechanism of their propagation in the econ- 
omy? Do cyclic oscillations form a constituent part of economic growth, or 
should they be regarded as deviations from a long-term trend? How should the 
government build an anti-cyclic policy, or, generally speaking is such a policy 
always necessary? What mathematical problems are encountered by researchers 
in their formalized description of oscillation processes in the economy? Exhaus- 
tive answers to these and many other important questions have not been found 
so far. 

From the point of view of synergetic economics, there are no economic evo- 
lutionary systems that always exhibit stability. An evolutionary system always 
experiences transformation effects of external and internal forces that are able 
to realize sudden structural changes, including cycles. We want to clarify this 
statement by the example of an analysis of behavioral properties of the phe- 
nomenon of a competitive interaction between different economic subjects ^5] . 

The notion of economic competition is a complicated and complex category: 
it is built on the basis of a variety of approaches. An overwhelming majority 
of researchers single out price competition and structural competition. In the 
first approach, the mechanism of competition is realized at the expense of price 
changes, whereas in the second approach, conditions of the production of goods 
undergo competition. Defining the essence of competition, one should take 
into account necessary constituents of the process that can be conditionally 
subdivided into three groups: behavioral, structural, and functional. 

Behavioral understanding of competition has been interpreted since the time 
of A. Smith as pair (without an agreement) competition for the most lucrative 
selling conditions that occurs between sellers (or buyers). At the same time, 
he considered price change to be the main method of competition. Later, the 
behavioral properties of competition improved in the direction of more accurate 
formulation of its objective and methods. 

The structural constituent of competition is formed by an analysis of the 
whole market or of its segments with the aim to determine the degrees of freedom 
of the seller and of the buyer. The functional filling of competition contains 
innovation, that is, a competition between the old and the new, etc. 

Before turning to a direct analysis of peculiarities of the models of economic 
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competition, it is necessary to point out the following: economic systems, as 
a rule, are far from their equilibrium, they are open to commodity and money 
flows, they have complex inhomogeneous structure and a regulation system of 
the endogenic medium under the influence of exogenic factors. Therefore, math- 
ematical formalization of the processes of economic competition is, in itself, of 
considerable diSiculty. This should be contrasted with the situation in physics, 
chemistry and biology where mathematics has already become a natural lan- 
guage of the description of the observed processes. With regard to speciflc prop- 
erties of economic phenomena, one refers specifically to mathematical models in 
economics. Here, by the model one should understand rather rude abstraction 
and idealization that represents mathematical formalization not of the develop- 
ing system of the economic space (market) itself but only of some qualitative 
and quantitative characteristics of the processes that flow therein. A general 
feature of many phenomenological models of the market economy is the pres- 
ence of autocatalytic (by analogy with biophysics and population dynamics) 
terms that determine the possibility of growth, of the facts of the appearance 
of unstable stationary states that facilitate the excitation of self-oscillation and 
quasistochastic regimes. 

Complex processes in the systems of market self-regulation are stipulated by 
the presence in the structural schemes of feedback contours (loops) (e.g., the 
"invisible hand" of Adam Smith); they are of both positive and negative types, 
which, in turn, predetermines the formulation of the problem of the study of 
structural stability of the considered objects and systems. In the equations of 
local competitive interactions, the feedbacks are described by nonlinear func- 
tionals whose character allows for the initiation of the excitation of complex 
dynamic regimes accompanied by corresponding attractors. The presence of 
nonlinear relations and of induced instabilities implies the use of the synergetic 
paradigm as a means of the description of competitive interactions in the eco- 
nomic medium. It is necessary to state that the classical linear principle of 
superposition loses its validity in a complex and nonlinear world represented 
by the market. In such a situation, it is impossible to argue that the whole is 
equal to the sum of its constituent parts. We should expect that the evolution 
of economic systems constitutes a specific transformation of all the participants 
of the market intcracti(m by means of the establishment of a coherent relation 
and of mutual adjustment of the parameters of their evolution. In this case, 
a nonlinear synthesis should be understood not as the unification of rigidly es- 
tablished, fixed objects, but rather as the unification of developing structures 
that are characterized by different economic "ages" and by "memory" and are 
positioned at different stages of the evolution. 

Thus, the complexity of the economic system is related to coherence. By 
coherence we shall understand the adjustment of the rates of business activity 
of the participants of the market by means of diffusive (mixing) dissipative 
processes that manifest themselves macroscopically as seeming economic chaos. 

The construction of a complex market competitive organization requires co- 
herent unification of the constituent substructures, the adjustment of time con- 
stants of their evolution. As a result of such a substantially nonlinear synthesis, 
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different structures get into a unified temporal world, that is, they acquire one 
and the same sharpening impetus and begin to function at one and the same 
economic rate. 

The development of the concept of the construction of the system of per- 
fect market competition leads to the understanding of the fact that arbitrary 
structures under arbitrary conditions or mutual relations, positioned at arbi- 
trary stages of development, cannot be included in a unified, complex social- 
economic formation. It seems that there exist a restricted choice of formations 
and of methods of the construction of a complex evolutionary whole. The se- 
lectivity and quantization of the methods of the unification of the parts into a 
whole is related to an imperative demand for existence in the whole temporal 
world. This is a natural, inherent basis for quantization in the process of inte- 
gration of complex dissipative economic systems. In the case when the economic 
subjects, integrated into a unified competition medium, differ in the sharpening 
impetus, they develop at different rates in the vicinity of a given singularity, 
which, in turn, provokes undesirable market imbalance. In the world economy, 
this means, for example, that the development rate, the standard of living, in- 
formation supply, etc. differ substantially for different countries and give rise to 
a dangerous difference in their potential. 

To restore the efficiency of market competition, it is necessary to observe 
certain topology of the "architecture" of the cross-relations. In other words, if 
the overlap region is small, the economic subjects will develop without "feeling" 
each other and will live in different temporal worlds. On the contrary, if the 
overlap is excessively large, the structures will quickly merge into strategic al- 
liances in the given market and possibly, may even form a unified dynamically 
growing structure with a growth limit equal to the market volume, which leads 
to degeneration of competition. 

Now, it is possible to proceed with the characterization of the mathematical 
essence of economic competition for several participants of the market. In this 
work, we have restricted ourselves to the case of the competition of two eco- 
nomic subjects in a single market. To construct kinetic equations of competitive 
interactions, it is necessary to make a number of assumptions that characterize 
the considered phenomenon at a qualitative level. It seems that the most im- 
portant issue is an analysis of the balance of the growth rates of the processes 
and factors that prevent positive development. Let us assume that the growth 
rate of each competitor depends on a potential increase in the volume of a cer- 
tain type of goods and on an unrealized possibility of growth for this type, as 
it has been the case for each isolated participant of the market in the absence 
of competition. However, the unused possibility of quantitative growth for this 
type of goods under the condition of mixing of flows of goods is a more complex 
quantity. It demonstrates the availability of a free place for this type of goods 
in the presence of goods expansion by another participant of the market. 

For two participants of market competition, we obtain a system of two au- 
tonomous nonlinear differential equations that resemble models of mathematical 
biophysics (of population dynamics) of the Volterra type [lOj, as well as their 
various modifications with various response functions [33j . 
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A theoretical analysis (including a bifurcation one) of such models in a suf- 
ficiently complete manner is given in the book by A. D. Bazykin 

Of special interest for the practice of economic forecasting is the formulation 
of criteria of closeness of the parameters of the system to dangerous boundaries, 
when, on crossing these boundaries, the system in a catastrophic way goes over 
to a qualitatively different state. In this case, the character of the dynamics of 
the volumes of goods in the market changes drastically: for example, a spas- 
modic transition from monotonous economic growth to relaxation oscillations 
takes place. Such boundaries of the regions of changes of the parameters of the 
considered dynamic system are called bifurcation boundaries. 

A special position in the studies of models of the competitive economy is 
occupied by processes that are characterized by cyclic behavior. The ascer- 
tainment of hidden periodicity, a search for the so-called " economic clocks" of 
various nature is always a valid issue in the studies of the problems of economic 
dynamics. 

Now, we shall discuss a mathematical model of economic competition that 
is described by a system of two ordinary differential equations with quadratic 
nonlinearity. Such systems appear in competitive dynamics when one uses a 
Taylor-expansion approximation of the second order for response functions in 
generalized Volterra models. The study of the issue of the excitation of self- 
oscillation regimes (limit cycles) is a very interesting and difficult problem of 
qualitative theory of differential equations. Up to now, there is still no solution 
to Hubert's sixteenth problem, posed in 1900: Find the maximum number of 
limit cycles and determine their mutual arrangement in a system of two differ- 
ential equations with quadratic nonlinearity. Among the main results for this 
system, we may note the following ^15j : 

1) a complete classification of their singular points (a node, a focus, a saddle, 
and a center) is given; 

2) a complete qualitative analysis of the systems with a center-type singular 
point is carried out. A topological classification of phase portraits is given, and 
a corresponding partition of the parameter space of such systems is performed; 

3) it is proved that limit cycles of quadratic systems are convex; 

4) limit cycles cannot surround a node- type singular point; 

5) a system that has an algebraic limit cycle in the form of an ellipse has no 
other limit cycles; 

6) a quadratic system that has a non-coarse focus and a phase straight line 
or two singular points with zero divergence has no limit cycles; 

7) a system with four singular points, two of which are focuses, with one of 
them being non-coarse, may have limit cycles only around one of the focuses; 

8) the maximum number of limit cycles generated by a focus or a center is 
equal to three; 

9) a quadratic system may have at least four limit cycles arranged as 3:1, 
i.e., three limit cycles around one focus and one limit cycle around the other 
focus; 

10) the total number of limit cycles in a quadratic system is finite. 
Three bifurcations of limit cycles are known: 
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1) the bifurcation of the generation (annihilation) of a limit cycle from a 
complex focus; 

2) the bifurcation of a separatrix cycle from a homoclinic or heteroclinic 

closed trajectory; 

3) the bifurcation of a multiple limit cycle. 

The first bifurcation is studied completely only for the case of quadratic 

systems: the number of limit cycles generated from the singular point is equal 
to three. For a system with cubic nonlinearity, the cyclicity of the singular point 
is equal to not less than eleven! 

As regards the second bifurcation, we may argue that, at present, there exists 
a complete classification of separatrix cycles, and the cyclicity of most of them 
is know. 

The third bifurcation is the most complicated one, and it is insufficiently 
explored. 

Unfortunately, all these bifurcations are of local character: when studying 

them, one considers only a certain sufRciently small neighborhood of a singular 
point, of a separatrix limit cycle or a multiple limit cycle and a corresponding 
sufficiently small neighborhood of the parameters of the system. 

A final solution of Hilbert's sixteenth problem requires a complete qualitative 
study of the system as a whole (both on the whole phase plane and in the whole 
parameter space), that is, global theory of bifurcations is needed. Besides, all 
local bifurcations of limit cycles should be joined together. 

Our detailed characterization of bifurcation properties of cyclic dynamics of 
the competitive interaction is by no means casual. Exactly here, the instabil- 
ity properties of the system with respect to small deviations of the parameters 
manifest themselves most strikingly. Only in nonlinear systems near bifurca- 
tion boundaries, qualitative differences in the character of the behavior of the 
considered object arc observed. An example of such restructuring of topology 
is provided by the transition from a stable aperiodic region to an unstable self- 
oscillation regime that occurs in a catastrophic manner. On the phase plane, 
this is illustrated by separatrices, i.e., lines that separate different attraction re- 
gions (attractors). In other words, if a small perturbation "throws" the system 
over a separatrix, it gets into the zone of influence of another attractor, which 
cardinally restructures the phase portrait. 

We have already studied in sufficient detail the qualitative peculiarities of 
market dynamics for two participants of competition. Certainly, this is only a 
particular case of complex organization of the economic system. It seems that 
the perfection of the market relations should be directed towards an increase in 
the quantity of the participants of the market. 

As is well known, the appearance in the market of a third participant of 
competition may initiate in the system a chaotic regime accompanied by the 
appearance of a new type of the attractor. This "strange" attractor radically 
changes the dynamics of the system of competitive relations, which substantially 
narrows the horizon of economic forecasting. Therefore, the most important idea 
that follows from synergetics is that stable development and the dynamically 
developing process of the evolution of the market necessitates a certain portion 
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of chaos, the spontaneity of development and self-organization, and a certain 
portion of external management on the part of state institutions that should be 
adjusted to each other. Both the two extremes, i.e., pure chaos, spontaneous 
market mechanisms of selection and the " survival of the strongest" , on the one 
hand, and total external management, full control and a protectionist policy 
vis-a-vis selected structures, on the other hand, are unacceptable. 
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